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length from left end of rotor to 
first major bearing station 

in. 
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acceleration . 
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— - 
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in. 
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total bearing force in x-dir. 
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radial clearance of journal 
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damping at major mass stations 
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mass station 

couple acting at i +h major 
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i n. 


xx 



NOMENCLATURE (Continued) 


SYMBOL 

DESCRIPTION 

UNITS 

ENX,ENY 

percent synchronous speed of 
harmonic force 

— 

ES 

equilibrium eccentricity 

In. 

ESU 

equilibrium eccentricity calculated 
from unbalance load 

in. 

EU 

unbalance eccentricity 

ml 1 s. 
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to support 
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relative shaft deflection at I™ 
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potential energy 

velocity 
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load parameter = '{Tlq /c/g 

support weight 
rotor weight 
journal weight 

bearing support weight 

critical speed of flexible rotor 
on rigid housing 

journal weight 

stability speed parameter, = <0/ /g/c 

total load on journal 

body fixed coordinate system 


UNIT 
rad, 
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lb. -in. 


in. 
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UNIT 
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rotor displacement (single mass 
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CHAPTER I 


ROTOR DYNAMICS: BACKGROUND AND STATEMENT OF THE PROBLEM 

l.l Introduction 

From the time the wheel was invented man has been plagued with 
problems of bearing lubrication and shaft dynamics. His ingenuity has 
produced the high-speed turbine rotors that power the giant jet liners 
of the present age and the vast array of process machinery that the 
American industrial community must keep in operation to meet the demands 
of our society. 

A typical example of a small high-speed rotating unit is shown in 
Fig. l.l. This is a cut-away of the Brayton-cycle turbine-generator 
unit that has a design speed of 50,000 RPM and is being developed for 
the NASA space program. The six-stage axial compressor and single stage 
drive turbine compose one integral shaft system with the generator rotor 
and drive turbine composing the second integral rotor-bearing system. 

Of great importance in such a unit is the design of the bearings and 
bearing housing mount structures. The present analysis will be concentrat- 
ing on methods of coup l.ing bearing «and rotor-shaft dynamics. 

A vast amount of information has been developed to date concerning 
the problems of lubrication or separately the kinematics and dynamics of 
machinery. The interest and means of solution for the combined problem 
of rotor-bearing system dynamics has been evolving since the early 1960’s, 
The earliest works in the fields of shaft vibrations and bearing whirl 
are important to the overall understanding of the scope of the problems 
of rotating equipment. The following background section will briefly 




review a portion of the earliest work but will concentrate on the 
present state of the art in the field of rotor-bearing dynamics, 

1.2 Background 

The study of rotor dynamics may be broken down into four distinct 
areas which have been treated in the literature. These are as follows: 

1) Critical speed analysis 

2) Unbalance response and the inverse problem, balancing 

3) Analysis of self-excited instabilities resulting from: 

a. Fluid-film bearings 

b. Internal friction damping 

c. Aerodynamic and hydraulic cross-coupling 

4) Time-transient simulation 

The well-known method of Holzer for finding the torsional natural 
frequencies was extended by Myklestad (I) for application to the analysis 
of aircraft wing structures and separately by Prohl (2) who formulated 
the problem specifically for flexible rotors. A tabular approach as 
advanced in their original works has been formulated for solution on 
digital computers and is currently being widely used for critical speed 
analysis of rotating shafts. 

One other major paper of the I940’s concerning critical speeds was 
the work of Green (3) which presented much data on the effect of gyroscopi 
on the single overhung disk, simply supported single disk and double disk 
rotors, and finally the infinite disk case. His work indicated the 
presence of both forward and backward critical speeds resulting from 
gyroscopic effects. 



More recent investigations (4-10) have made slight extensions of 
the excellent work of Green and Myklestad-Prohl to include the effects 
of shear deformation and specialized cases such as the two-bearing 
machine with an overhung disk. Yamamoto(ll) in 1954 presented a major 
contribution to the study of critical speeds in which he presented both 
theoretical and experimental observations on the phenomena of forward 
and backward whirl due to gyroscopic effects. Most current critical 
speed codes calculate only the synchronous critical speeds of forward 
precession. 

Closely related to the problem of critical speeds is the analysis 
of unbalance response. A major contribution to steady-state response 
was the method of Lund (12,13) which uses a modified Myklestad-Prohl 
technique to solve for unbalance response. A detailed discussion of this 
approach is presented in Section 2.4. The analysis can include known 
speed dependent linear fluid-film bearing characteristics and support 
flexibility and damping. Analysis of multi-mass flexible rotors has also 
been conducted by Tang and Trumpler (14) and Koenig (15). The effect of 
support damping is very important from the viewpoint of reducing vibration 
amplitudes but even more important for the reduction of forces transmitted 
through the bearing and to the support (16,17,18). 

The rotor modeled in Fig. 1.2 has been studied (19) as an extension 
of the early work of Jeffcott (20) and others (2i,22,23) to include 
support flexibility and damping. The critical speeds of the resultant 
two-degree of freedom system are shown in Fig. 1.3 where the critical 
speeds divided by the rigid support critical are plotted versus the 
ratio of support to shaft stiffness, K, for various values of the ratio 
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FIG. 1.2 SCHEMATIC DIAGRAM OF A SINGLE MASS ROTOR ON DAMPED 
ELASTIC SUPPORTS 
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CRITICAL FREQUENCY VS SPRING RATE 
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of support to rotor mass, M. 

The unbalance response of this simple model can be used to illustrate 
several important aspects of rotor-bearing response characteristics. For 
example the absolute rotor response is shown in Fig. 1.4 (a) for various 
values of the ratio of support to rotor damping. A lightly damped support, 
c = 0.01, produces two distinct resonance zones in the given response 
range. The accompanying phase shift of the deflection with respect to the 
unbalance is shown in Fig. 1.4 (b) for the response of 1.4 (a). The 
lightly damped system is observed to go through a 180° phase shift at the 
first critical and another phase reversal as the second critical is 
passed. This phase reversal is suppressed for higher support damping as 
indicated in Fig. 1.4 (b). The force transmitted to the bearings is 
shown in Fig, 1.4 (c) where the addition of proper damping clearly would 
reduce the level of the forces transmitted for the bearing support. The 
phase shift to 180° and continuing up to 270° at the second critical is 
a very important concept to note and can be of great importance in the 
analysis of experimental data. The damped support is clearly advantageous 
for reduction of the forces transmitted in or near the critical speed 
regions of operation. 

The inverse problem of unbalance response is the ever present problem 
of balancing of rotating equipment. The method of Thearle (24) is still 
being used for field balance procedure. Two basic approaches to the 
problem of balancing flexible rotors have been presented in the literature. 
The first, modal balancing, consists of balancing at a speed corresponding 
to the first mode, then successive trials at each of the higher modes to 
reach the desired speed of operation. The other method known as the 
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PHASE ANGLE (PHI 2) [DEGREES] AMPLITUDE RATIO (X2/E) [DIM. 3 



FREQUENCY RATIO (W/WC) 


(a) ABSOLUTE ROTOR MOTION WITH A TUNED SUPPORT SYSTEM FOR VARIOUS VALUES 
OF SUPPORT DAMPING, K=M=1,A=10 



FREQUENCY RATIO (N/WC) 


(b) PHASE ANGLE OF ABSOLUTE ROTOR MOTION RELATIVE TO UNBALANCE FOR VARIOUS 
VALUES OF SUPPORT DAMPING 

FIG, 1.4 AMPLITUDE, PHASE, AND FORCES TRANSMITTED FOR A TYPICAL 
SINGLE MASS FLEXIBLE ROTOR ON DAMPED SUPPORTS 
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influence coefficient approach utilizes information obtained from added 
trial weights and response data at each of several runs to produce a 
matrix of coefficients which is used to predict the resultant unbalance 
in the rotor. 

A new approach presented by Little (25) uses the influence coefficient 
technique coupled with the application of linear programming to arrive at 
the optimum response by proper choice of correction weight size at limited 
locations along the rotor span. The problems associated with balancing 
will not be directly treated in this presentation but the literature on 
this subject is a valuable portion of the theory of rotor dynamics (14, 
26-34 ) . 

The literature contains many references to the observation and 
probable cause of self-excited instabilities. The earlier reports by 
Newkirk (35), Robertson (36), and others (37,38) on internal friction and 
hydrodynamic bearing instabi I ity are well referenced by Gunter (39), 

Ruhl (40) and also in (41,42,43). A review paper by Newkirk in 1957 (44) 
points out one distinguishing feature between ’'cramped-shaft whirl" (ie., 
internal friction damping) and oil-film whirl. The former may become 
resonant above the critical speed whereas the later does not become 
resonant at speeds below twice critical speed. In an earlier presenta- 
tion (45) Newkirk had distinguished between certain phenomena that may 
occur as a result of journal bearing instability. 

(1) A "half frequency whirl" may occur at any running speed 
the whirl frequency being slightly less than half the 
running speed 

(2) A "resonant whirl" which occurs when the running speed 

is equal to or greater than twice the natural frequency of 
the rotor for transverse vibration at that speed. 



the oil film 

(4) A subharmonic vibration of an unbalanced rotor, the 
frequency being a definite fraction of the running 
speed. 

i 

Morrison and Paterson (45) presented an approach to describe the stability 
threshold of symmetrical two-bearing flexible rotors by one equation and 
thereby reducing the criterion of other investigators such as Hori (47), 
Holmes (48) and Pinkus and Stern I i cht (49) to the general form 


w o to 9 
(_V < P (Jl) 2 - Q 
0) — 0) T 
9 


D.i] 


where 

id = velocity of rotor 

u> n = natural frequency of simply supported rotor 
P,Q = functions of the running condition 
and 

,V 2 _k me _ c 

a) m W " 6 

9 

where 

c = radial clearance of bearing 

6 = static deflection of rotor centerline under gravity 

loading 

Ruhl (50) has recently presented an analysis of the stability of the 
short journal bearing having shaft flexibility and plots “/“g versus 
eccentricity for various values, of shaft static deflections (6/c). The 
introduction of the flexible shaft reduced the stability boundary below 
the value for rigid rotor analysis. (See Fig. 6.6). 
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A recent work by Choudhury (43) has examined the effect of linear 
support flexibility and damping on the plain journal bearing. Numerous 
stability maps were presented for various ranges of support mass and 
flexibility. The analysis then investigated the stability of a multi- 
mass rotor system by producing the characteristic polynomial of the 
system by numerical methods. The stability study of a particular rotor 
system was presented which investigated the effect of support flexibility 
on the system and results verified that an optimum damping exists for a 
given support stiffness to improve the stability of the rotor system. 

Two text books on the subject of rotor dynamics have been published 
which give detailed derivations of the equations of motion of the 
fundamental rotor models. DImentberg, in 1961 published a text (51) 
which treats flexural vibrations in rotating rotors and considers various 
aspects of critical speeds, balancing, stability of rotors operating 
near critical speeds, and the influence of shaft internal friction and 
shaft asymmetry on stability. Experimental results are reported in the 
text which add to its importance in the rotor dynamics literature. The 
text of Tondl (52) investigates instability in rotor systems due to 
internal friction damping, unequal stiffness, and oi I film journal 
bearings. This extensive treatment is also supplemented by reports of 
experimental investigations carried out by the author. Tondl has one- 
fourth of the total volume dealing with the effect of journal bearings 
on stability. Experimental results for plain cylindrical, elliptical, 
lobed, grooved, and flexible element bearing are presented and discussed. 

Tondl points out one major fact that should be emphasized. For 
machines of high power rating, measurements made directly on the rotor 
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are the only reliable guide for the correct assessment of the danger 
of vibrations. Pedestal mounted pickups are not reliable for determining 
the level of vibration on the rotor shaft. 

Another mechanism creating instabilities in turbomachinery has been 
reported by Alford (53). Aerodynamic forces arising from labyrinth seals 
and blade-tip clearance in axial compressors and turbines can cause a 
resultant driving force component and should not be neglected in 
simulation studies of real machinery. The analysis of Choudhury 
discussed earlier in this section had examined this form of instability 
in the example rotor system presented in that work. 

Kerr (54) reported experimental observation of a flexible, damped, 
air lubricated journal bearing in which a region of whirl could be 
passed through to a stable operation zone. Gunter (39) had earlier 
presented a thorough treatment of the effect of internal friction on the 
flexible support rotor system. Support symmetry and absence of support 
damping was shown to lower the stability threshold of the rotor system. 

An added observation pointed out the fact that the stability threshold 
predicted from linearized theory using for example the method of Routh, 
did not necessarily represent the limit of safe operation. The growth 
rate of a given instability may be very low and operation may be possible 
above theoretical stability boundaries. This fact points out the 

t 

necessity for the transient solution technique which can examine the 
nonlinear growth rate of an instability predicted by a linearized 
stabi I ity analysis. 

Analog computer simulation of transient whirl for a symmetrical 
rigid rotor in cavitated short journal bearings was reported by Jennings 
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and Ocvirk (55) as early as 1962. Steady-state orbits for balanced rotors 
were presented which had whirl rates of approximately one-half but no 
orbits which encircled the origin were produced in that analysis resulting 
from limitations of the analog system. Gunter (39) also presented analog 
computer produced orbits for both linear bearing characteristics and 
nonlinear fluid-film bearings. Orbits showing combined effects of 
synchronous unbalance response and half-frequency whirl instability were 
presented. The nonlinear case formed finite limit cycles whereas the 
linear systems were unbounded. 

Reddi and Trumpler (56) in 1962 examined the stability of the 360° 
full journal and the 180° partial-film bearing. End leakage factors 
were applied to the film force expressions and the resulting equations 
of motion linearized and examined for stability about the equilibrium 
position by the Routh criteria. The complete equations of motion were 
programmed on a digital computer and the resulting orbits presented for 
the 360° journal. Reddi also demonstrated that if cavitation is not 
included the bearing will be unstable at all speeds. 

Caste 1 1 i and McCabe (57) presented transient simulation of a rigid 
rotor supported in tilting pad gas bearings. Results of two case studies 
are given with experimental verification and excellent agreement. 

Elrod, McCabe, and Chu (73) reported an approach to determine 
stability boundaries of gas-bearing systems by a combination of the 
transient orbit and linearized equation methods. The system linearized 
characteristics are calculated using gas-lubrication theory and these 
linear characteristics are then' used in transient simulation to deter- 
mine stability in the small. Digital computer transient simulation was 
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also reported by Badgley and Booker (58) who presented stability 
boundaries for the rigid, balanced rotor system having cavitated journal 
bearings and several transient whirl orbits illustrating the balanced 
system response. Kirk and Gunter (18) used the cavitated short bearing 
theory developed in a fixed coordinate system to study the effect of 
unbalance and external loading on the stability of the horizontal rotor 
system. The values of the instantaneous forces transmitted were given 
for selected cases in addition to absolute whirl rates. The maximum force 
transmitted was indicated for all whirl orbits in the presentation. The 
stabil izing effect of certain levels of unbalance on the otherwise un- 
stable vertical journal was also investigated by the time-transient 
digital computer simulation program. 

Akers, Michaelson, and Cameron (60) reported an analysis of finite 
length journal bearings including friction effects and out of balance 
loading. Stability boundaries were compared to the work of Badgley and 
Booker (58) and was found to be in fair agreement. Shapiro and Colsher 
(61) presented an analysis of a resiliantly mounted hybrid journal 
bearing and an elastic, gimbal mounted Rayl eigh-step thrust bearing. 

Both time-transient and step-jump dynamic analysis was used for the 
results on the stability and response of the systems studied. The 
step-jump approach is used to calculate bearing system characteristics 
at a given position. Then these linear equations are placed in 
characteristics to diminish the time that would be required for a com- 
plete nonlinear transient solution and is similar to the work of Elrod 
reported earlier (73). 
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The development of pseudo-transient techniques (73) for the predic- 
tion of journal bearing performance in internal combustion engines has 
allowed accurate determination of minimum film thickness, bearing loads, 
and overall system performance. The stability and transient motion of 
a vertical finite three-lobed bearing with cavitation was investigated 
by Falkenhagen, Gunter, and Schuller (62) and compared to experimental 
work and the analytic work of Lund on the three-lobe bearing. Transient 
orbits for both balanced and unbalanced systems were presented using an 
approximate finite bearing solution which represented a time saving 
factor of 100 over the complete numerical finite-difference approach to 
the solution of the governing Reynolds' equation. 

The time-transient approach is not limited to the study of lubrica- 
tion problems and the power of the approach has been real ized in many 
sciences. The transient analysis of multi-link mechanisms has been 
implemented (63,64) as well as the transient solution of the dynamic 
behavior of a closed-cycle gas turbine system (65). The numerical 
solution to any time dependent system of equations has become a usable 
tool for the engineer-scientist due to the speed of the present generation 
of digital computers. 

Shen (66) has presented a formulation for the flexible rotor using 
influence coefficients and the assumptions of small displacements and 
a torsional ly rigid system. Provision was made in the derivation for 
analytical representation of nonlinear bearings. The derivation also 
assumed that the bearing mass was negligible at the two major bearing 
stations and a chord approximation was used to define the gyroscopic 
torques and thus reduce the number of parameters in the solution. 
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However, at the same time additional fictitious stations must be specified 
for the desired accuracy of angular deflections. The procedure outlined 
by Shen also involves a transformation to rotating coordinates which is 
not required or thought to be advantageous in the following analysis. 

Childs (67) has also presented a formal derivation of the dynamics 
of flexible rotating system in terms of a spinning reference coordinate 
system using vector mechanics. These equations are then reduced to the 
fixed cartesian coordinates by a transformation. This formulation takes 
into account bearing constraints which may be nonlinear, but no provision 
has been made for bearing support flexiblity. The approach by Childs 
requires that the system eigenvalues and eigenvectors first be obtained 
to allow a modal coordinate transformation and omission of modes higher 
than running speed if desired. 

These formulations for the transient solution of flexible rotor- 
bearing systems are excellent approaches to a very difficult problem. 

The following analysis presents a derivation of the equations of motion 
by yet another approach which readily allows the inclusion of nonlinear 
time dependent bearing, seal or external shaft forces. Fluid-film 
bearings with the option of flexible, damped bearing support structures 
is included in the simulation model . Linearized gyroscopics are included 
in the derivation and the effects of acceleration on the transient 
response is included in the analysis and is an option in the computer 
code developed from the theory. The following analysis investigates 
steady-state and transient response of both rigid and flexible rotors 
and stability programs for special rotor systems are also presented and 
used in connection with the transient analysis. 
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1.3 Statement of the Problem 

The purpose of this analysis Is to present the formulation of the 
time transient solution of a flexible rotor on nonlinear fluid-film 
bearings that may be used for design studies of real rotor-bear i ng 
systems including the determination of bearing support characteristics 
for optimum system performance. 
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CHAPTER II 


FLEXIBLE ROTOR DYNAMICS 

2.1 Description of the Simulation Model 

The large class of rotating machinery operating near or above their 
first bending critical must be analyzed under the general field of study 
known as flexible rotor dynamics. A typical rotor in general has multiple 
bearings and the rotor shaft may have disks or impellers located inboard 
and outboard to the main bearings. 

The cross-section of the complex Brayton-cycle turbine-compressor 
as shown in Fig. 2.1 is typical of the high speed units which can be 
studied for dynamic behavior by analytical simulation. This particular 
unit has a six stage axial compressor supported between the main bearings 
with the single-stage drive turbine located on the right overhang. A 
thrust pad is located on the front of the compressor as noted in the 
diagram. 

The continuous rotor shaft may be regarded as a series of concentrat- 
ed mass stations connected by massless elastic shafts. For practical 
purposes the gyroscopic effects should be included in the analysis of 
overhung rotor stations but may be neglected in most cases for rotor 
stations inboard of the main bearings. Although the notation flexible 
rotor is appropriate, the deflections considered are small in comparison 
to the dimensions of the rotor and simplified equations for the rotor 
shaft mechanics may be incorporated into the analysis. 

Fig. 2.2 represents an idealized rotor system reduced to the model 
used for the analytical description. This analysis will consider the 
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rotor to be supported on two main bearings which may in turn be elastically 
mounted. Any number of additional bearings or seals may be incorporated in 
a given analysis and is included as an option in the computer code (see 
Appendix B). 

The rotor station deflections will be denoted by the deflections 
relative to the main bearing centerline (u j ,v. ) and also by the correspon- 
ding absolute deflections (Xj,y. ). These are shown in Fig. 2.2 for the 
y-coordinate deflection in the y-z plane. The relative and absolute 
angular displacements for the i^ rotor station are also shown in Fig. 2.2. 
Angular displacements of each disk will be denoted by the angles 0^, 
which is a rotation in the y-z plane about the negative x-axis and 0 

X 

which is a rotation in the x-z plane about the positive x-axis. 

A cross-section of the ? mass station is given in Fig. 2.3 which 

indicates a possible whirl configuration for the eccentric disk at 

some instant of time. The absolute and relative deflection nomenclature 

for the station geometric center, o^, is il lustrated and the displacements 

of the bearing support and journal are denoted as (X , T ) and 

s i s i 

(X.,Yj) respectively. These deflections, which represent the rigid rotor 
centerline and support centerline at this instant in time, are referred 
from the first bearing location as follows: 






- a 


•) x [> t 


a] 


C2.il! 


V 


X - X 
S2 Si 

( ) xl*| - a] 


b - a 


[ 2 . 2 ] 
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Similar expressions hold for the y -coord i nate deflections referred to the 
i +h station. 

2.2 Equations of Motion 

The equations of motion of multi-mass flexible rotor are derived in 
the following sections. The equations are developed in fixed cartesian 
coordinates and include the effects of gyroscopics and acceleration. 

The relative deflections of the mass stations are assumed to be small so 
that linearized gyroscopic moments can be used in the equations. The 
flexibility influence coefficients for the rotor shaft are assumed to 
remain constant for the order of approximation being considered. 

2.2.1 Gyroscopic Moments 

Consider an isolated mass station with symmetry about the z-axis as 

shown in Fig. 2.4 with the polar and transverse moments of inertia given 

by 1^ and 1^ respectively. The standard Eulerian angles are shown in 

Fig. 2.4 and are denoted as 0 the nutation angle of the shaft axis 

referenced from the fixed Z-axis, <j> the precession angje as projected in 

the X-Y plane and referenced from the X-axis, and the local body 

rotation about the z-axis. Also shown on the figure are the angles 

(6 ,0 ) which are the projections of the nutation angle 0 in the fixed 
x y 

X-Z and Y-Z planes respectively. 

For a rigid body the angular velocity may be expressed in terms of 
the Eulerian angles and referred to the body-fixed coordinates (x,y,z). 
The expression for this relation is readily available in most dynamics 
texts but a rigorous derivation is usually not included. A derivation 
of this expression using unit vectors and orthogonal transformations is 
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FIG. 2.4 TYPICAL GYROSCOPIC STATION SHOWING EULER IAN ANGLES 
AND THE ANGULAR DEFLECTIONS © x AND 6 y 
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presented In Appendix A. 3. If the unit vectors n x> Py> and n z lie along 
the (x,y,z) coordinates, then the angular velocity may be expressed as 
fol lows: 


= (S sin ip - <j> sin 6 cos iJO 

t (e cos i|> + $ sin 0 sin ip) riy + ($ cos 0 + ^) ri 2 C2.3] 

The equations for the inertia torques are most easily derived by 

energy considerations and the application of the Langrangian. 


± ,2k-) . lk_ = o 

dt V 3q| ' 


where 


L = T - V 

T = kinetic energy 

V = potential energy 

Q = genera I i zed forces 

q = generalized coordinates 


C2.4] 


The kinetic energy is expressed as the sum of the translational 
kinetic energy and that due to a rigid body rotation of S$, That is 


T> = * «l *J S ij 


T 2 • 


i I,j «,j 


.(no cross products of inertia) 
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Concentrating on the contribution due to rotations it is desired to obtain 

the equations of motion in terms of the angles 0,0 and angular velocity 

x y 

0 ). 


T 2 = i I T (Q*r}) 2 t -J- I (n • n ) 2 + I • n ) 2 
z T 'x T y p z 


or substituting Eq. 2.3 for ft 


T 2 = ■§• Iy(0 2 + $ sin 2 0) + £ Ip (ij) + cos 0) 2 [2.5] 


This expression is in terms of the Eulerian angles and must be 


transformed to the 0 , 

A 

for small 0 (from Eq. 


0^ coordinates. 
2.3). 


The following relations are valid 


( 1 )^ — <J> t \f) [2.6a] 


FoC infinitesimal rotations the angles may be treated as vectors and 
the following expressions may be written. 


0^ = 0 cos <J> 

= 0 cos <t> -t 0$ sin <(> 
0^ = 0 sin <J> 

0 ^ “ § sub 4 > t 04 » cos <J> 


[2.6b] 

[2.6c] 

[2.6d] 

[2.6e] 
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For small 9 the kinetic energy T 2 reduces to 


T 2 ~ % Ij(9 2 t (j»9 2 ) t {j- Ip(^) + <M I ~ 2 ~ + ...)) 2 C2.7] 

As shown by Yamamoto (II), 

9 2 t 0 2 = 0 2 t 0 2 $ 2 C2.8a] 

x y 

Vx ' ®x 9 y = £2. 8b] 

Hence, 


T 2 = £ I T (0 2 + ®w 2) + ^ I + ♦ " 4>9 2 / 2 3 2 

1 x y p 

T 2 = 1 + ®y 2) + £ *pE w2 " a> ^ 82 + 4> 04 /4H 

or 

T z = i Iy(9 x 2 + 9 y 2 ) + i I p Co) 2 + O)(0 x 0 y - 0 y 0 x )D C2.9a] 

For a torsional stiffness K 

v = i K(e 2 + 0 2 ) [2.9b] 

x y 


Application of Eq. 2.3 to 2.9 gives the desired equations for the 
generalized moments M x , M y where the rotations (-0 y ,0 x > correspond to 
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2.2.2 Equations of Motion - Rotor Stations 

The equations of motion of the flexible rotor may be obtained by 
application of the flexibility influence coefficients. A discussion of 
the method of obtaining these coefficients for a stepped round rotor 
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shaft is presented in Appendix A.I. The influence coefficients are 
defined as 


= the relative deflection at the i station due to 

_l~h 

a unit unidirectional load applied at the j station. 

= the relative deflection at the i** 1 station due to a 

"t”h 

unit torque at the j n station. 

i"h 

= the relative angular rotation at the I station due 

to a unit unidirectional load at the j +h station. 

= the relative angular rotation at the i^* 1 station due 

to a unit torque at the station. 


The relative deflection of the I major mass station can be expressed 
in terms of the flexibility influence coefficients as follows: 


u i “ij P Xj + 6 ij C y. 


[ 2.153 


V i °ij P Yj " 6 U °xj 


[ 2.163 


6 xj " *'IJ P Xj + Y ij C yj 


[ 2 . 17 ] 


0 = d». . P - y. . C 

Yj 'J Yi 'J x, 


[ 2.183 


where the resultant loads and couples are composed of both inertia 

J J 

loading and effective externa.1 loads created by unbalance, internal 
damping* cross-coupled aerodynamic loading, and any reactions from 
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fluid-film bearings or seals at the particular station. 

The sign convention for the loads and couples is shown in Fig. 2.5 
which gives a top view of the i^ mass station. Also shown are the 
absolute and relative translational and angular deflection notations. 

The inertia loading acting at the i +fl mass station is readily derived 
by differentiating the position vector, R. (see Fig. 2.3). 

-> -> ->■ “*■ _ __ 

R. = x.i + y.j + e. cos(wt+(j>. ) j + e. sin(a>t+^.)j C2. I9] 

R- = -jp = (x,y). + (-eu sin(utt<f>), eu cos(ut +<j>) ) . C2.20!] 

^ •• •• * 

Rj = (x,y).+ (~e u 2 cos (cot +6 ) -e u sin(utt<j>), 

-eu 2 sin (utt<|> ) + eu cos (ut+4> ) ) . C-2 . 2 1 U 

The inertia loading including unbalance is then given by D'Alembert 
principle as 

inertia) _ _ m . x (acceleration). = (-mx,~ my). + 

+ (meu 2 cos(ut+4>) + meu sin(ut+<j>), 

meu 2 sln(utt<j>) - meu cos(utt<|> ) ) . H2.22J 

The expressions for the gyroscopic inertia torque and angular 

stiffness were previously derived CEqs. 2.13-2.14) and in terms of 

the absolute angular rotations (© x ) these may be written as 

J J 

f o 1 1 ows : 
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FIG. 2.5 FORCE AND MOMENT SIGN CONVENTION AT 1 th ROTOR 
STATION 



C - ^' ner *^ a ^ 
Xj " x 


. - ( I_0 — o)I 6 - I w0 + KQ ) . 

J Ty px py yj 


[2.23] 


C = (^inertia) 

y j" y 


. = (-l_0 - wipe - i I u0 - Ke ). [2.24] 

j T x K y ■ p x x j 


Combining this expression with the representation of internal damping (68), 
cross-coupling (53,69), and allowing for other nonlinear forces that 

A A 

might occur described in general terms by F x (x,y,x,y) and F^(x,y,x,y), 

the components of P. may be expressed as follows: 

si 


P = (-mx - cx - kx t mew 2 cos(wtt<j>) t 

X j 

A 

t mew sin(wtt(j>) - Clu - wClv - Qy + F (x,y,x,y)) 

x j 


[2.25] 


P = (-my - cy - ky t mew 2 sin(wtt<j>) 

y j 

A 

- mew cos(wtt<f>) - Cl v + wClu + Qx - W + F (x,y,x,y)). [2.26] 

y j 


Rewritting Eqs. (2.15-2.18) and using the notation 


X j 


X j 


= -m.x. t F 
J J x. 

[2.27] 

= -m.y. + F 
J Y J Yj 

[2.28] 

= I T 0 + H 

T j y j X j 

[2.29] 

= — I T 0 t H 

T J X j y J 

[2.30] 
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then 




» U 


x, U 


V V- mv + Vj - VVWj 


+ Vj 

[2.313 

' Vj 

[2.323 

+ Vj 

[2.333 

• H ) . 

x J 

[2.343 


The equations of motion of the rotor may be expressed in matrix 
form for convenience. On the following page is the matrix formulation 
for two major mass stations. 

The formulation for larger systems is simplified by defining the 
following variables. 


a ll rn l a l2 m 2 
0 * 21^1 ^ 22^2 


«m m r 


5 A 


(n) 


a nl m l 


a m 
nn n 


3llI T 1 &l2l T 2 ‘ * * ^in^T 


^T, 


= B 


tn) 


Bni* 


Ti 


e i T 

nn * 
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Then the equations may be expressed very simply as 


A (n) 0 (n ) Q (n) g (n) 
Q (n) A (n) B (n) Q (n) 
Q (n) c (n) D (n ) . Q (n) 
c (n) Q (n) Q (n) Q (n) 


where 


.(n) 


“I 

>>■ 


’p(n)" 




X 


y Cn ) 

S 

p(n) 

y 


e£ n) 


^(n) 


y 


y 


e<"> 


^(n) 


L X J 


X 


: (n) _ 


lj F x . 
j 

+ 3. .H 

,J y j 

.F 

+ e .h 

nj X j 

nj y. 

lj F y. 

“ p ij W x. 

J 

, .F 

- e .h 

nj y j 

nj Xj 


C2.35] 


•©* 


,(n) 


lyi + +>j F y. ’ Y ij H x. 
• J J 


■0 + A .F - y .H 

y n Y nj yj ’nj Xj 


& 


(n) _ 


S + ♦lJ F x. + *lJ H y 

• J 


•0 + <j> .F + y .H 

xh nj Xj T nj y. 


This formulation presents the equations of motion which are in 
the desired form for application of standard integration procedures. 
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The method of solution employed in the computer code MODELJ will be 
discussed in Section 2.3. 

If the system experiences an acceleration u, then one additional 
relation must be given to fix the instantaneous angular velocity. That 
is 

. do) 

“ = dt 

Hence 

dio = didt 


or 

t 

a). = a), + / wdt 

0 t. 

For a constant acceleration rate then 


u> + = w + t uj (t - t Q ) C2.36D 

o 

The angular velocity can thus be calculated for each time step of 
the integration process. The external torques necessary to supply 
this constant acceleration is given by Eq. 2.12. 

2.2.3 Equations of Major Bearing Support Stations 

The previous equations are all that are necessary when considering 
a flexible rotor on two rigid bearing supports. If the bearing or 
supports are allowed to have flexibility, then additional information 
must be obtained to solve the system. Taking the moments of a I I 
external and inertia forces and torques acting on the shaft about the 
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first bearing station (see Fig. 2.6) : 




n 

(b-a) + l C 
i = l y i 



[ 2.373 


Hence, the total resultant- force component acting on the shaft at the 
second major bearing in the x-di recti on is given by: 


F 

x 2 


n n 

- I c - l P x U, - a) 
i = l y i i=l x i 1 


(b - a) 


[2.383 


This resultant force is composed of the bearing reaction denoted as 
(B.F.) and the inertia load at the second major bearing. That is 


F 

X2 


-m. x. + (B.F.) 

J 2 *2 x 2 


[2.393 


Therefore, solving for the journal acceleration x, at the second major 

J 2 

bearing station yields 


m, 


{(B.F.) 


+ JT~ — -r- 

X 2 (b - a) 


n 

1 I C y 

i = i V l 


n 

t l P. U, - a))} 


i=i 


'i 


[2.403 


The bearing reactions can also Include unbalance excitation acting 
at the journal. In a similar fashion the acceleration for the first 
bearing may be expressed as 





m, 


{(B.F.) 


*1 



n 


Me 

t=i 


n 


+ I P (b - V>>> 
i=l x i 1 


ko 


[2.413 




FIG. 2.6 REPRESENTATION FOR BEARING FORCE AND MOMENT 
BALANCE 


k } 



The acceleration of the bearings in the y-coord i nate are likewise found 
to be given by the following expressions. 


n n 


*J 2 - — {<B - F - ) y 2 + E-hr ( -,I,°*, tl.Vh - am 


i=l I i=l ' i 


and 


n n 


'Ji • ±- + btT <+ T, C x. + T, P y. ,b - V» C2.43] 


1=1 I 1=1 'I 


If the bearings are supported by flexibly mounted structures, the 
equations of the supports are readily deduced using the representation 

f 

given in Fig. 2.7. A force balance on the support gives 


If 


““ m 3 
support s s 


= -<B.F.) + (S.F. ) 


C2-443 


where 

(S.F.) = support force 

a g = absolute acceleration of the support In any given coordinate 
direction. 

The bearing and support forces have been left in the general terms 
since the bearing type assumed dictates how the reaction forces are 
expressed. These forces may be either linear or nonlinear functions 
of displacement, velocity, and time. For a more detailed analysis of 
the bearing and support forces refer to Appendix B where the types in 
the computer code MODELJ are discussed. More detailed discussion of 
fluid film bearings is presented in Section 4.2 and discussion of the 
fluid film squeeze damper is presented in Section 6.3. 
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SUPPORT SYSTEM FORCE BALANCE 


* 3 ' 


I 


2.3 Method of Solution 


The equations derived in Section 2.2 can be solved for the time 
transient solution once the inverse of the modified flexibility matrix 
is obtained. Expressing the equations of motion given in Eq. 2.35 as 
follows for simplicity: . 

CCXA] DO = 00 C2.45II 

where 

DG ~ column vector of absolute acceleration of the 
major mass stations. 

OG = column vector of relative displacement as 
represented in Eq. C2.35U 

CCXA] = 4n x 4n matrix given in Eq. C2.35] 

Then, premultiplying both sides of Eq. 2.45 by the inverse of the 
4n x 4n matrix CXA gives the absolute translational and angular 
acceleration rates which is expressed as follows: 

ex: = ccxat'lu: [2 . 46 ] 

With the acceleration values known at time t, the bearing accelera- 
tions can be obtained from Eqs. 2.25 and 2.26. As for any time integratior 
of second order equations, the displacements and velocities at time t 
must be specified for all the coordinates being solved. The actual method 
of solution may be chosen from several available formulations. For the 
type solutions given by unbalance response the modified Euler method has 



been used successfully (57,59,62) and Is especially desirable when It Is 

necessary to calculate nonlinear hydrodynamic bearing forces at each step 

of the solution. A modified Euler method which is equivalent to a second 

order Runge-Kutta procedure can be utilized at no extra time loss in the 

integration of the bearing forces. A 4 +h order Runge-Kutta procedure 

can be used but becomes time consuming since four calls of the bearing 

forces are required for each step forward in time of the solution. The 

numerous predictor-corrector integration schemes such as the Milne, 

Adams-Bashforth, and Hammings (18) are even more time consuming and have 

"t"h 

a greater tendency to exhibit numerical instability. The 4 n order 
Runge-Kutta procedure is used in the time transient programs listed in 
Appendix E and G. 

The most important information needed to produce readily recognized 
trends Is the starting conditions. For the linear model the gravity 
loading can be removed from the solution and any external loading can 
be investigated from the static equi I ibrium position. For nonlinear 
solutions (bearing and support forces) the gravity loading cannot be 
ignored and the Initial conditions given may produce transient motion 
that the system would never encounter in actual practice. The following 
section is presented to review work that has been completed on the 
steady state response of flexible rotors and a computer code modified 
from (12) is presented in Appendix H that may be used to obtain initial 
conditions for the transient solution program. 

2.4 Steady State Response 

The transient solution as proposed in Section 2.3 is indeed very 
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Time consuming and hence very costly to produce even with the tremendous 
speed of the present generation of digital computers. Poor or misjudged 
starting conditions for the time transient solution may produce misleading 
results since the system would have to encounter a transition period to 
attain the correct phase relation for the many variables being integrated 
simultaneously. 

The equations of motion as presented in Section 2.3 can be solved 
for the steady-state for the case of assumed linear bearing characteristics. 
By assuming solutions of the form 

q. = A. cos(oit) t Bj sin(wt) [2.473 

the solution requires that a system of equations be solved. For the 
assumption of N point-mass rotor stations then a 4N x 4N matrix must be 
solved. However, if the shaft properties are isotropic then only a 
2N x 2N system need be solved. The introduction of support flexibility 
increases the system of equations to a 8N x 8N matrix representation 
and increases the solution time per case by approximately a factor of 
64 for a simple 2 station rotor. The inclusion of gyroscopic moments 
at each rotor station in addition ;fo support flexibility increases the 
system to a I2N x I2N representation and is 216 times as difficult as 
a two station point-mass rotor with isotropic constraints. Recall that 
the transient solution requires only a 4(N-2) x 4(N-2). matrix inverse 
be taken even when gyroscopics are included in the solution. For point- 
mass assumption the transient solution requires only that an (N - 2) x 
(N - 2) matrix be inverted. The (N - 2) expression represents the total 
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number of shaft stations in the simulation minus the two major bearing 
stations. A more efficient procedure for obtaining the steady state 
solution has been reported by Lund (12) which uses a modified Myklestad- 
Prohl technique to obtain the steady state solution. 

Without reproducing the entire analysis of Lund the basic equations 
for this type analysis will be sumarized in the following discussion. 
With reference to Fig. 2.8 the moments and deflection relations can be 
derived by applying simple beam theory for the shaft sections and noting 
that 


M . . = M» + L V 
nt| n n n 


[2.483 


d 2 x M 


By integrating the equation from simple beam theory, g— j and using 

the above expression for the moment 


9 n + , ’ + C / L " & H ^ + C /" #3 »„ 

z=o O 


[2.493 


n+l 


dz-i r n zdz 

e r" 


= X + L 0 + L [ / - / 

n n n n ' EI J { 

*?n 7/J7 L n 2 j y 

+ V L n / I r-/ tt-> 


[2.503 


For the case of constant shaft properties (El constant for each different 
shaft, section), then 


M . , * M' t L V 
nt| n n n 


[2.513 
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i»8 





Z2.521 


e 


x 


n-H 


n+I 


6 n + ir M n + ^r V n 


L* 

n 


L 3 

n 


x + L 0 4* — T M + £ ct m 
n n n 2Ei n 6Ei 


[2.53] 


The moment change across a mass station includes the gyroscopic 
action and angular spring rates of the bearing. The shear change across 
a mass station includes the effect of inertia loading, bearing reactions, 
and unbalance loading. Detailed equations for these contributions are 
given by Lund (12). Harmonic motion is assumed and the bending moment, 
shear, angular and lateral deflections may be expressed as 


= M cosut + 14 si nut 
xc xs 

[2.54] 

= V cosut + V si nut 
xc xs 

[2.5511 

- 0 cosut + 0 si nut 
c s 

[2.563 

= x cosut + x si nut 
c s 

[2.573 


and for the y-direction the variables M y , V y , <f>, y are used in similar 
equations. The resulting equations needed to describe the rotor 
reactions include the change of moment across a mass station and shaft 
length (8 equations), change of shear across a mass station (4 equations), 
and the angular and lateral displacement along a shaft length (8 equations) 
for a total of 20 equations to step across one mass and shaft section. 

For a beam with free ends- the bending moment and shear are zero at 
the extremeties, or they can be calculated for an assumed deflection and 
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given support characteristics. The unknown quantities are the angular 
and lateral displacement components at the left or first rotor station. 
Using the superposition principle, each unknown is given a unit value 
with all the others set to zero and the twenty equations as mentioned 
above are used to calculate the bending moment and shear at the right 
rotor end. An additional calculation is made applying the given un- 
balance loading with the other unknowns set to zero. The result of 
the superposition of the applied unit loads produce a set of equations 
which may be solved for the unknowns at the first rotor station. This 
procedure is best described as follows: 

Let 9 c i = I, set 0 sl = <j> cl = <J>si = x c l = x sl = Yci = Ysl 

= u^ = u = 0 (u^ n and u^ n are unbalance loads) 
Calculate the residuals at the right rotor end. Call them 


M» . 


M* . ; M* 

xcy,l 

1 ’ xsy, 1 

r ycr,r ysr. 

V 

• V 


xcr, 1 

I* xsr,l 

1' ycr,r ysr. 


Let 8 = I and set the other variables to zero. Calculate 

si 

the residuals and denote them as 


1’ 


M’ 

m’ 

xcr, 2' 

xsr,2' 

ycr ,2 

ysr, 2 

f 

V • 

V 

V 

xcr, 2’ 

xsr,2' 

ycr, 2’ 

ysr, 2 


Repeat the procedure for, the eight displacements and then 
apply the unbalance loading and call the residuals 
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1' 

M' 

M* n : 

M' 

xcr,9' 

xsr,9 

ycr ,9 

ysr ,9 

, 

V n ; 

V ol 

V 

xcr,9' 

xsr,9' 

ycr, 9’ 

ysr, 9 


This produces a set of influence coefficients which allows the 
unknowns to be solved. The solution can be expressed as 
fol lows: 





9 

9 

9 

9 

9 

9 

9 

9 


[2.583 


The unknowns at the first station can thus be determined and the 
response all along the rotor can then be calculated for a given speed. 

Lund's original analysis included an iterative approach to the 
solution of the gyroscopic co.nTri.but Ion but the linearized gyroscopic 
equations as presented in Section 2.2.1 have replaced this iterative 
approach in the computer code given in Appendix H. One additional 
modification has assumed circular synchronous precession which restricts 
the bearing and support characteristics to be symmetrical. This allows 
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a tremendous savings in solution time and rotors having unsymmetrica I 
bearing and support characteristics may be studied by this program if 
the gyroscopic coupling is small in comparison to the other shaft loads 
and reactions. The results of this computer program, LUNDJR, may be 
plotted by a companion program LUNDJRP which is given in Appendix H 
also. A sample case is given in the following rotor system. 


Example 2. I 

Total rotor weight 

Center span approximate weight 

Journal weights (each) 

Diameter of shaft sections 
Length of each shaft section 
Youngs modulus 
Symmetric bearing stiffness 
Bearing damping 

(Rigid bearing support pedestals assumed) 
Unbalance at center station 


1301 lb. 
676.0 lb. 
312.5 lb. 

6.0 in. 

34.5 in. 

30 x 10 6 lb/ in 2 
448,000 lb/ in 
501 lb-sec/ in 


5.4 oz.-in. 


This system may be visualized by letting K^ x = = °° in Fig. 5.1 

and Kj x = Kjy = 448,000 Ib/in while Ci x - Cjy = 501 Ib-sec/in. The 
equivalent shaft stiffness is 280,000 Ib/in for this example case study. 

Fig. 2.9 represents the undamped critical speed mode plots 
corresponding to the rotor of example 2.1. It is seen that the rotor 
first critical speed at 3300 RPM is symmetric with the largest deflec- 
tion at the rotor midspan. This should be the predominant mode 
excited by unbalance placed at the rotor mid-span or center. The 
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THREE- MASS 
IFFNESS= 




second critical speed is a conical rigid body mode and will not be 
excited by a centrally located unbalance. The largest amplitude occurs 
at the bearings and hence should be a highly damped mode for typical 
bearing damping as given in Ex. 2.1. 

The third critical speed shows considerable shaft bending with the 
largest amplitude occuring at the bearings. Note that the amplitude 
at the bearings are out of phase to the rotor amplitude at the midspan 
section. This mode may be excited at the rotor center if the bearing 
damping is not excessive. The steady state response of the center mass 
station is shown In Fig. 2.10 and the corresponding phase angle with 
respect to the unbalance Is given in the joining plot. In the speed 
range shown one critical speed is evident at 3300 RPM and the phase 
shift through 90° to 180° is shown in the adjoining plot. The bearing 
amplitude and phase is shown in Fig. 2.11 where the phase angle is 
observed to cl imb beyond 180° to more than 300°. The second critical 
is not excited by the given unbalance and the third critical is com- 
pletely damped out of the response curves. 

It should be quite evident that this computer code alone could 
be used to optimize a rotor bearing support system over a given speed 
range. Once a given configuration was decided upon the transient 
response analysis could be used to determine the sensitivity of the 
design to shock loading and to determine the stability or instability 
of the system. 

The fourth critical speed shown in Fig. 2.9 is a direct consequence 
of including the gyroscoplcs of the lumped rotor mass stations. A point 
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ROTOR- BEARING SYSTEM 



mass model calculation would only have produced the first three critical s 
shown In Fig. 2.9. The steady-state response given In Figs. 2.10 and 
2.11 did not cover the speed range of the fourth critical speed and 
unbalance at the center station would not have excited the mode even if 
the range was extended to beyond 14,600 RPM. This mode could be 
excited by an unbalance moment such as that caused by a skewed disk. 

The inclusion of gyroscopics in general will lower the critical speeds 
if the transverse moment of inertia Is larger than the polar moment of 
inertia. However, In most applications to compressor and turbine disk 
locations the polar moment of inertia is larger and the inclusion of 
gyroscopics will tend to stiffen the shaft and hence raise the critical 
speeds. 
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CHAPTER III 


DYNAMICS OF A TWO STATION ROTOR HAVING 
ONE SECTION OVERHUNG 

3.1 Description of the Simulation Model 

One class of industrial rotors may be represented by a model as 
suggested in Fig. 3.1. (a). This class of machines has a massive mid- 
span section and an overhung section supported in fluid film or rolling 
element bearings. If the effect of the bearings can be neglected and 
gyroscopics are neglected on the midspan, then the model is reduced to 
the system as shown in Fig. 3.1 . ( b ) . This representation has a point 

mass at the midspan station and a massive disk at the overhung section. 

3.2 Equations of Motion 

The system being analyzed must be represented by six equations of 
motion. The midspan characteristics are denoted by the subscript I 
while the overhung station will be denoted by the subscript 2. The 
equations for the deflections are given in Chapter II and may be written 
as fol lows: 


x 2 x 2 


'xj + °'2 p x 2 + 6 l 2 °y 2 

D.a 

Vi + “ 12 P y 2 ' 6 l 2 °x 2 

C3.2] 

J ' + ClooP + 3?9^ 

X} 22 x y 2 

2 

[3.3] 

5 + ao^P 3? 2^ 

yi 22 Y2 22 x 2 

[3.4] 

3 + <booP + Y 99 C 

Xi *2 Vl 

[3.5] 
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where 


p Xl =-mix 1 - cxSci - kjxj - CljX! - uCliVx - Q iyi 

- k x, - k yi - c xi - c yi 
xx i xy 7i xx 1 xy 71 

p yi * ~ m lYl ~ c lYl " k lYl ~ CI lYl + wCliXi + QiXj 


" k Vi “ k Xi - c Vi - c Xi 

yy/l y x 1 yy / 1 y x i 


P •— "^ 2^2 C 2 X 2 ^ 2^2 ^^ 2 X 2 *“ CoCl 2 Y 2 “ Q 2 Y 2 

x 2 

p „ = -rn 2y2 - c 2 y 2 - k 2 y 2 - ci 2 y 2 + u£l 2 x 2 + Q 2 x 2 


C 

*2 


I T 0 
T y 2 


' mI P 8 x 2 


+ KoQ 


Y2 


G 

Y2 



“ I P 9 V2 



[3.7] 


[3.8] 

C3.9] 

[3.10] 

D.ll] 

[3.12] 


This formulation requires that the flexibility influence coefficients 
be known and the resulting equations are coupled in the acceleration 
terms as were the equations in Chapter II. The method and computer code 
used to Investigate the stability of these equations will be explained 
in the following section and the results of the time transient solution 
is presented in Section 3.4. 

3.3 Stability Analysis 

The equations of motion given by Eqs. 3. 1-3.6 may be examined for 
stability by applying the Routh criterion (69). This procedure requires 



that the system characteristic equation be known. If solutions of the 
following form are assumed: 


x l = Aje*^ 

» . At 

Vl = a 2 « 

a At 

x 2 = A 3 e 

A *+ 

y 2 = A 4 e 


0 x 2 = ^ 




the equations of ^motion may be readily written as a six by six matrix 
formulation. Substituting the above expressions into Eqs. 3. 1-3.6 and 
collecting terms produces the matrix formulation presented on the 
fol lowing page. 

A computer program, OCSTB, was written to expand this matrix which 
has elements quadratic jn A to obtain the characteristic equation, 

j 

examine the resulting characteristic equation for stability by the 
Routh criterion, and to fi$d all the roots of the resulting twelth order 
polynomial (69). The real part of the roots indicate the stability of 
the system (positive for instability) and the imaginary part gives the 
natural frequencies of the system, A listing and explanation of the 
input required for OCSTB are presented in Appendix C. 

The stability analysis produces the natural frequencies of the 
system. As an example consider the following rotor system. 
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Example 3. I 

Wi = I .86 lb. 

W : 2 = 1.86 lb. 

I p = 0.00542 lb. -In. -sec . 2 
I T = 0.00274 lb. -in. -sec . 2 
overhung span = 3 in. 
bearing span = 12.75 in. 
shaft diameter = 0.375 in. 

E = 30 x I0 6 lb. /in . 2 

The influence coefficients were calculated from the procedure 
described in Appendix A. The results are as follows: 

a u = 1.483 x 10 " 3 a 12 = -1.047 x 10 " 3 812 - 3.489 x I 0~ 4 

a 21 =-1.047 x 1 0 ' 3 a 22 = 1.622 x 10" 3 8 22 =-5.923 x I 0 ~ 4 

<|> 21 = 3.489 x I0 -4 <fr 22 = -5.923 x I0 “ 4 y 22 k 2.49 x I0 -4 

The natural frequencies were calculated and are given in Fig. 3.2. 
The solid lines are the forward critical speeds and the dashed lines 
give the backward criticals, (3,43). Note that the scale on the ordinate 
is broken and the gyroscopic critical is very high indeed. The 
gyroscopics increase the two lower forward criticals as the speed in- 
creases. The synchronous forward critical speeds are approximately 
2750 RPM and 6250 RPM as taken from Fig. 3.2 The critical speeds and 
mode shapes of synchronous forward precession were also obtained by a 
program (70) using the Myklestad-Prohl technique. These results are 
given in Fig. 3.3-4 where it Is noted that the bearing characteristics 
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CRITICAL SPEED ft, RPM x 10' 


IO 


FIG. 3.2 CRITICAL SPEEDS OF FORWARD AND BACKWARD PRECESSION 
FOR A TWO MASS ROTOR SYSTEM 
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10000.00000 

FIG. 3.3 MODE SHAPE OF AN OVERHUNG ROTOR AT THE FIRST FORWARD 
CRITICAL (OVERHANG ■ 3IN. ) 



10000.00000 

FIG. 3.4 MODE SHAPE OF AN OVERHUNG ROTOR AT THE SECOND 
FORWARD CRITICAL (OVERHANG * 3 IN.) 





were chosen to approximate rigid supports. The critical speeds were 
given as 2745 RPM and 6284 RPM by this computer code and it was con- 
cluded that the two programs were in agreement. The mode shapes given 
could be made smoother by simply taking more points along the rotor 
shaft. The rotor of Example 3.1 was modified slightly and will be the 
basis for the analysis of the remainder of this chapter. 

Example 3.2 

This rotor is the same as the rotor of Example 3.1 with the 
following exception: 

overhung station = 5 in. 

The influence coefficients for the rotor are: 

an = 1.483 x lO -3 a 12 - -1.744 x I0~ 3 0 12 = 3.489 x I0 -3 

a 2 i = -1.744 x JO -3 ot 22 = 5.079 x I0“ 3 8 22 = -1.159 x 10" 3 

<\> 2 i = 3.489 x I0‘ 4 <J> 22 = -1.159 x I0" 3 Y22 = 3.176 x I0" 4 

The resulting critical speed plot is given in Fig. 3.5 where it 
is observed that the extended overhung section has lowered the two 
forward bending criticals N T , N TT to 1825 RPM and 4960 RPM respectively. 

if ii f 

The corresponding mode shapes are given in Figs. 3.6 and 3.7 Also 

indicated on Fig. 3.5 Is the effect of internal friction damping acting 

at station one and at station two. With a value for internal damping 

of 0.2 lb. -sec. /in. acting at the rotor midspan, the system threshold 

occurs at a shaft speed of 2500 RPM and whirls at 75% of running speed 

which is observed to be the lowest critical N, at this operating speed. 

L f 

For the same internal damping at the overhung station the system 
threshold occurs at 1890 RPM and whirls at 9556 of running speed which 
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CRITICAL SPEED H.RPM x |$ s 



FIG. 3.5 CRITICAL SPEEDS OF FORWARD AND BACKWARD PRECESSION 
SHOWING INFLUENCE OF INTERNAL DAMPING ON WHIRL RATIO 
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FIG. 3.6 MODE SHAPE OF AN OVERHUNG ROTOR AT THE FIRST 
FORWARD CRITICAL (OVERHANG * 5IN.) 



4966 RPM 


0 

CM 
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FIG. 3.7 MODE SHAPE OF AN OVERHUNG ROTOR AT THE SECOND 
FORWARD CRITICAL (OVERHANG = 5IN.) 


corresponds once again to the lowest forward critical of the system. 

At higher speeds the whirl remains at that frequency corresponding to 
the lowest critical speed. The results for a rotor speed of 3000 RPM 
produces a whirl which is 62$ of the running speed as shown in the 
computer results given in Table 3.1. 

The rotor specifications are printed at the top of the table with 
the corresponding influence coefficients. The matrix for the system 
is shown and corresponds to the elements of the system as shown on 
page 60. Each element of the 6 x 6 matrix has three coefficients with 
the first one being the coefficient of the quadratic term, X 2 , the 
second the coefficient of X and finally the constant term. The roots 
to the resulting characteristic equation are also shown with the real 
and imaginary parts, then the imaginary part in RPM and finally the 
fraction of running speed for each solution. 

3.4 Transient Whirl Analysis 

The simple model of Example 3.2 was chosen to test the transient 
program which was formulated in Chapter 2, Of particular interest is 
the effect of the inclusion of gyro scop ics on the results of the time- 
transient analysis. The rotor shaft specifications used in the analysis 
are presented in Table 3.2 along with the approximated rotor specifica- 
tions which correspond to the rotor of Example 3.2. A gyroscopic station 
on the overhang and a dummy seal (bearing) location are also specified 
in Table 3.2. For the purpose of simplifying the initial conditions 
required a zero gravity field is assumed. 

For a rotor speed of 3000 RPM the results of the stability program 
with no internal friction damping are as given in Table 3.3, These 
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TABLE 3.1 SYSTEM CHARACTERISTICS OF A TWO MASS ROTOR 
CONFIGURATION HAVING INTERNAL DAMPING AT THE 
OVERHANG 
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(UNITS ARE INCH-PUUNO-SECCNOI 


KOTOR KID SPAM SPECIF! CAT X OKS 

hi i,46 csi .« *oi 

KSi « a. co cn * a. so 

oi_t a^ia lejljs cuuii 

kxx * o.ea cxr * o.oo 

KXT *_™. 0.03 CYX » 0.90 

kyx • a. 09 ctr * a.aa 

KTX-* a. oa 


ROTOR OVERHUNG SPECS 

N2_* 1.96. 

* XJ» * .0QS420 CS2 * .61 

XU _*OOZT40 ..Cl 2. * JI .0Q. 

K 52 * 0.00 

«U 5 JUJLG 


INFLUENCE COEFFICIENTS 

All-* 1.483D-C3 U2 * -1.74fcD-03_ 512 _»..3.A*M-8* i 

A21 * -1.7*60-03 AZ2 * 5.0190-03 bZt * -1.1590-03 

pel t ezlj? mz * ?«!?&, or, fl* 


MATRIX FOR OVERHUNG CO**PR 

I 7.1451CD-06 1 0, “l -8. 483733^0 eT l“~ 0. I 9.S5946J-C7 I 

I. 1.402800-05 I 0. I -1*244000-05 1 8. __ I 0. I 

t l.g<iqo30*qo 10 . i a. i o. I o. x 

-ah r — — rr: — r: — 1 ~~ --^ x : 

X 0. I 7.145100-06 I tf. I -4.493TT0-J« I 6. I 

I 1.442800-05 I 0. ... I. ._-l*744O3O-05 X -5.94J940-44 1 

10. I 1. 18.0 000+ 08 I 3. I 0. X C. X 

-I—- “I . .--I -I—-- ---i.-- ... I- 

I -8.-03730-06 I 3, X 2.447393-G5 I 3. \ X -3.175663-00 I 

j: -i.. 7450 CD- 0 5 I 0. I — 5.379530-C‘? I C. I 3. L_ 

0. X 0. I 1*303003 *i*J X 0. X C. X 

10. I -8.403730-06 I 0. l' 2.4-7390-35 I 6. I 

I 8. .1 _-l. 744000- pS X 0. X__„ _ 5.07903 J-j5 I 1.973440-33 X 

X 0. 10. I 8. I 1. v06t»33+3Q I 0. X 

" 5ls=a 1 - — s~=ssl — r--— i-- -i- 

I 1.681230-06 X 3. X -5 .S6-42D-46 X U. X 4.7022-J-V7 X 

.1. 3.489000-06 I .0. I - 1. 159033-03 I . 0. X . 3. I 

I 0. X 8. X 6. X 0. X 1.003030*03 X 

— X — ——I — ----- — — I— rr I — — — — x — — I- 

10. I 1.681230-06 I 0. I -5.544820-36 I w. I 

J. fi* I 3.469000- S6 I 0. I -1.1590 30-35 1 -5. 53/930-34 X 

10. I 0. 16. 10. I 0. X 


.0. I 

5.94CO3J-0* A 
t. * t 
i ■ 

9.554360-0? I 
j . > 

C. i 


>i._ 973 - 4 J - tf 3 : 
•J.l76uo j-66 . 


5 ! 4 C? 930-84 I 

6 . I 


• . * l 

1. Suk.l) j* Dl 1 


C HARACTERISTIC JETTON 

GIVEN ll* FORM AO * Al+L ♦ A 2 *L * # 2 ♦ ... # 15 n**U» 

l.ToMli+00 ” 1.312360-04" ""bT 527620^05 5^74 5t 0-09 “ 1.285970"- 09 474 23740- 1^'~7. 369*0 J-~l5 

1*074330-19 1.366440-20 . 2.704370-26. _3.3t6630r2L_— ^ 7.2C240-33 2.972310-3*.^ , 

CQEEFICXEHTS I JL 3 £Jt£SS£ NP^AUIEP 

i.OOQCQD+OO 0 . 301000+00 1 . 600440+07 _ . 1 . 305 Q 1 D+ 3 S 6 . 593400+13 5 . 184230*14 _ 3 . 556130+19 _ 

2 . 134690+20 6 . 205510+24 2 . 456970+25 3 . 149930+29 6 . 332860+29 4 . 425540+33 

SRPN ■ 3000.01 KONTBI * 0 . MOUNT « „0 . ■ ... _ 

CXI* 0.00 CI 2 * 0.00 , • 


THE ROOTS OF THE CHAR. £0. ARE AS F0U06S _ 


-9.686A767D-01 
-9. 64687670-01 
-1. 05264550*30 

-1.79502600+02 
1. 79502600+02 
. , -1^9465 0173+02 _ 

-1714,1 

1714.1 

-1858.7 

RPM. _ 

KPN. 

-.5714 

.5714 

-.6196 

-1.05264550*00 

1.9464 U 1 70 + 02 

1458.7 

RP*f . 

.619o 

-1.01763520*00 

-5,10379650+02 

-4873.4 

PPM. 

-1.6246 - _ _ _ . . 

-1.01763520*00 

5.10379650*02 

4473.6 

RPN. 

1.6244 

-1.03275990*00 

-5.17761730*52 

_ -4944.3 

KPH, 

. -1.6481 _ . _ 

-1.03275990*00 

5.17761740+02 

4944.3 

RPH. 

1*0461 

-S.gflflCEMQ-g? - 

-2. 4597 76fl J*03 

-23489,2 

OPH, 

„ -7,6297 _ ... - 

-5.20800610-02 

2 * 469 7 76 00* 33 

23469.2 

KPN. 

7.4297 

-2.67292950-02 

-3.C5869630+03 

. . _ -25208*4 

RPT. 

-9.7361 

-2.67292950-02 

3.95469690+03 

29208.4 

APN . 

9.7361 


TABLE 3.3 SYSTEM CHARACTERISTICS OF A TWO MASS ROTOR 
CONFIGURATION SHOWING STABLE RESPONSE FOR THE 
CASE OF NO INTERNAL DAMPING 


lb 








results predict a stable operating condition. Since the gyroscopics 
are not predominant at this low speed a point mass rotor simulation was 
run using the code MODELJ . The rotor specifications for this assumption 
are given in Table 3.5 (see page 76). 

The Initial conditions were given to excite the first mode response 
with a synchronous initial whirl rate. Fig. 3.8(a) represents the 
transient response of the overhung disk for five cycles of running speed. 
The motion has an average whirl of nearly 60 % and is stable as predicted 
by the stability analysis. The midsection response given in Fig. 3.8(b) 
indicates the presence of a higher mode developed in the x-di recti on. 

The running speed is between the first and second critical s and a 
combination of the modes is to be expected for free response of the 
system. 

If internal damping is introduced at the midsection the stability 
analysis indicates an unstable whirl of 62 % as given in Table 3.4. The 
rate of growth of the instability should not be too great due to the 
relatively small real part of the root. The specifications for the 

1 i 

transient simulation are given in Table 3.5. The transient orbits for 
synchronous first mode excitation are given in Fig. 3.9. The first 
five cycles for the overhung and midsection are given in Fig. 3.9(a) 
and 3.9(b) respectively. The orbits indicate an approximate whirl of 
60 % with definite indication that the motion is orbiting out into an 
unstable pattern. The continuation for cycles six through ten are 
shown in Figs. 3.9(c) and 3.9(d) with definite indication of an 
< unstable system. The 60 % whirl has predominated the response at the 

'' ■■ . \ '■ .j. 

midsection as indicated in Fig. 3.9(d). 
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FIG. 3.8(a) TRANSIENT MOTION AT OVERHUNG STATION 
GYROSCOPIC MOMENTS ( H=0 . 05 ) 

(b) TRANSIENT MOTION AT MIDSPAN STATION 
GYROSCOPIC MOMENTS (H=0.05) 
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1 . 0 b 

0.00 
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.01 
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(MX « 

0.00 

CAY • 

0.00 

—RXV « 
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CTA • 
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K VS » 
*W 1 » 

0.06 

8.18 

CJf • 

6.86 

ROTOR 0 VERMUM 6 SPECS 



1.06 


_ 

t* • 

•065620 

CS 2 * 

*11 

at.*^ 

.602760 
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- . 8 * 08 . - 

«S 2 • 

. 07 A- 
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INFLUENCE COEFFICIENTS 

—Ail » 1.4030-33 Al2 * -1.74*0-03 612 m 3.4090-04 
021 • -1.7440-03 A22 * 5. 0790-03 422 « -1.1540-03 
4-WQgQ a B22 .. K Ti . tsar i T Aj .^y.M ,<.W60-fl%— 


WAT* IX FOR OVcRNUNt COHPR 

I 7.145100-0* X 0. • 2 -4.403733-06 t 0. X *.»59A4Cr-fl7 1 

X J. 113B4D-G* 10. X -1.744000-05 l 0*. X 0. X 

X 1.000080*00 1 9.314730-02 X 0. X 0. X 0. X 

Jiflr -.-.r--.— r.7-- I -.T.r-r r-: — r----— - It-- r- — - I- - r -r r- - X- — — ■ 

10. 1 7.145100-0* X 0. 2 -8.403730-06 I 0. X 

X 0. X 3.11 3880-06 X 0. . I *1.766000-05 X -5.9*0480-0* X 

X -9.316730-02 I 1.000000*06 X 0* I 0. X ft. X 

X *0.463730-96 x 0. I 2.4472 W-05 X 6. X -3.175663-0* X 

fl, ,-L 3/3.600 r.OS-1 fl. J 0. X_ 

X 0. I -1.695790-01 X 1.800000*00 X 0. X 8. X 

— X— — — — — — -X— — — — — — — -X— 1 — — -I— -—I- 

X 0. X -0.463 730- 66 X 0. X 2.447390-05 X 0. X 

X 0. X -3.662400-06 X 0. .1 5.079300-05 X 1.97340^-03 I 

X 1,095790-01 I 0. X 0. X 1.000060*00 X O, X 

^lr^r*-r,-=.r. — It — — — rr.X-: 

X 1.601230-0* I 0. X -5. 506*20-06 X 0. X 8. 7 <2 22*0-07 X 

X 7.326900-05 I 0. I *>1.159603-05 X 0. X 0. 2 

I 0. X 2.192210-02 X 0. X 9. X 1.600000*00 I 

— x— X— ————.< 1- — •— I— — — — — — — X-— — — — 1— 

X «. X 1.601230-06 I 0. X -5.506020-36 X 0. I 

— ,1 U - I 7.324900-05 I 8. r -1.154000-35 Z_ _-S.« fl/93P-IK I 

X •2.192210-02 X 0. X 0. I 0. X 0. X 

— X—^———— -***-X— —r I- — — - X — —I — — X- 


0. X 

5. 5606*0-06 X 
6. X 

'17555*6 J-0> l 
0. X 

0. X 

0. X 

.Zl. 9736.80- 03 I. . . .. 
0. X 

-3.175060-06 X 
0. I 

C. I 

ir ■ x 

5. **07930-8* I 
0. X 

... — 

0.702260-07 i 
UL* _J 


.000600*60 X 
— — — — X— 


_CH4.<%^reiu>ix.c_u«.ATj:p»i , • 

feXVtN IH FOKrl AO ♦ Al*U ♦ A 2 *l **2 ♦ .. . Al 3 *t** 12 » 


1.060*00*00 7.505290-04 6. 560650-05 ' ' 3.346 500-08~ 1.292170-09 " 4.077070-13 7.397100-15 

1.176570-10 1.367160-20 2.901090-25 3.317160-27 .1.0922*0-32 _ 2. 072 410-3* 

XQMTX6UTO_iN„m^ — — 

1.660000*00 9.131100*01 1.660710*07 1.430560*04 6.597200*13 5.677580-15 3.469500*19 

1.9*7610* 21 6.235*60*26 1.613900*26 3.169730*29 3.623660*30 4.0*7630*33 

SRPM ■ 3600.01 ROMTai * 0 KOUNI • 0 _ .. _ . 

cn * .20 CX2 * 0.00 


JME ROOTS OF THE 

CHAR. eq. are as fouoks 



-8.050*1210*00 

-1.00082870*02 

-1719.7 RP?I. 

-.5732 

-8.05461210*00 

1.00402870*92 

1719.7 RPH. 

• 5732 

8.2 6 02*920-01 

-1.4*696500*02 

-1659.2 RPM. 

- —6197 

8.2602*920-61 

1.96*9*6wti*Q2 

1059.2 RPM. 

.6197 

-7.97929360*00 

-5.17*942 40*02 

-49*1.7 RPH. 

-1.6672 

-7.9/929J0O*OO 

5.17*90290*02 

4941.7 RPH, 

1.6472 

•3.0271/310*01 

-3.02717310*01 

-5.09594270*02 

-6066.3 RPH. 

-1.6221 

5.69593270*02 

4666.3 RPH. 

1.6221 

_j^5 t 223jl«*O^02 

-3* 06469690*04 

-29208.4 RPH. 

-9.7361 

-5.223014*0-02 

3.05869660*03 

29200.4 RPH. 

9.7361 

-1.14072930-01 

-2.46377**0*03 

-23689.1 RPH. 

-7.8297 

-1.10072940-01 

2.46977*60*04 

23689.1 RPH. 

7.8297 


TABLE 3.4 CHARACTERISTICS OF TOO MASS ROTOR SYSTEM 

WITH LIGHT INTERNAL DAMPING AT MIDSPAN (CIl= 
0. 2LB-SEC/IN) 
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FIG. 3.9(a) TRANSIENT MOTION OF OVERHUNG STATION WITH INTERNAL DAMPING 
AT MIDSPAN STATION 

(b) TRANSIENT MOTION OF MIDSPAN STATION FOR FIVE CYCLES OF MOTION 





SEGHEMT I N-INIT. IN-FIHRL I T-FINRL 
(BPH> <RPH) I 
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FIG. 3.9(c) TRANSIENT MOTION OF OVERHANG STATION FOR CYCLES 
FIVE THROUGH TEN (CIl =0.2) 

(d) TRANSIENT RESPONSE OF MIDSPAN FOR CYCLES FIVE 
THROUGH TEN (CIl =0.2) 





The inclusion of the gyroscopic properties of the overhung section 
should not drastically alter the transient response from that obtained 
for the point mass simulation. With the same time step as for the 
previous solution but including the gyroscopics at the overhang the 
transient response indicated from the program is shown in Fig. 3.10. 

The initial response is smooth but small cusps are developed which grow 
into inner loops and finally the response is dominated by an approximate 
1 0/1 whirl rate. The cusps and inner loops are actual solutions of a 
simple gyro (71) but the large jumps in the later portion of the 
solution are clearly inaccurate. There are no energy producing mechanisms 
In the system and hence the physical considerations do not allow this 
type response (72). 

The system- has been excited by a frequency ten times running 
speed and Fig. 3.5 shows that the gyroscopic critical Is very nearly 
ten times running speed for the case of N = 3000 RF*M. The initial 
conditions on the angular displacements of the overhung were specified 
as zero and hence induced small excitations of the order of 1 0/1 of 
running speed. The time step of the solution was not small enough to 
accurately track these excitations and they eventually became dominant 
due to numerical difficulties. The solution at the end of the third 
cycle of running speed was used for initial conditions with the step 
size reduced to H - 0.01 (radians). This gives 628 steps per cycle or 
for a 1 0/1 whirl there would be 63 steps per whirl revolution. A 
sustained whirl of 1 0/1 should result from such starting conditions and 
the solution should be bounded. The results for one cycle of running 
speed is shown in Fig. 3.11 where the expected sustained 10/ I whirl is 
evident. 
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GYROSCOPICS SHOWING INSTABILITY DUE TO RELATIVELY 
LARGE TIME STEP (H = 0.05) 

(b) MIDSPAN RESPONSE SHOWING EFFECT OF NUMERICAL 
INSTABILITY AT OVERHANG DUE TO GYROSCOPIC NUMERICAL 
INSTABILITY 
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ftlSIQ BCAPINO SUPPORT SYSTEM 

TABLE 3.7 SPECIFICATIONS FOR ROTOR OF CASE NO. 112.7201 
(GYROSCOPICS INCLUDED) 


SEGMENT N-INIT. N-FINH. T-FINffl. 
C BPH) <RPM) <B30) 


u_ 



a 



(S1IM) UI(U 
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FIG. 3.11(a) OVERHANG RESPONSE CONTINUE FROM CYCLE THREE 
OF CASE 112.7201 WITH TIME STEP REDUCED TO H=0.01 
(b) MIDSPAN RESPONSE WITH H=0.01 FROM CYCLE 
THREE OF CASE 112.7201 SHOWING STABILITY OF 
NUMERICAL SOLUTION 






The next case considered was to start the solution with the initial 
conditions as in Fig. 3.10 but using a reduced time step of H = 0.01. 

The rotor specifications are indicated in Table 3.8 and the results for 
two cycles are shown in Fig. 3.12. The overhung station is observed to 
have a slight excitation of approximately 1 0/1 as indicated by the whirl 
path from cycle I to the end of this solution orbit. The reduced time 
step has not allowed the excitation to grow as it had in the previous 
case with the larger time step. 

The next example case considers the influence of unbalance on the 
resulting transient including gyroscopics. The unbalance specifications 
are as given in Table 3.9 which indicates 0.002 oz-in. unbalance at each 
of the two major mass station and 180° out of phase. The response for 
five cycles is given in Fig. 3.13 where the gyroscopics are observed 
to once again excite the solution and alter the orbits much the same 
as before. The step size in this case was the larger value of H = 0.05. 
A smooth orbit is obtained up to the end of the third cycle of motion. 
Additional translational damping at both stations was introduced as 
indicated in Table 3.10 with the results obtained shown in Fig. 3.14. 

The initial response is suppressed but the gyroscopic influence persists 
to enter the solution in the fourth revolution. The point mass model 
as indicated in Table 3.11 produces the response shown in Fig. 3.15. 

The overhung station has a very regular symmetric pattern and would 
approximately repeat itself for the next five cycles as would the mid- 
section. 

The next consideration was to again reduce the time step and start 
the solution at the end of cycle three for the case of Fig. 3.13. 
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SPECIFICATIONS (CASE NO. 1217.7101) 
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FIG. 3.12(a) RESPONSE OF OVERHUNG STATION SHOWING NUMERICALLY 
STABLE SOLUTION INCLUDING GYROSCOPICS (H=0.01) 

(b) RESPONSE OF MIDSPAN STATION (H^O.Ol) 
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TABLE 3.9 ROTOR SPECIFICATIONS FOR UNBALANCE RESPONSE 
AND LIGHT DAMPING 



SEGMENT I N-INIT. N-FINfiL T-FINflL 
(BPh) <RPM) (fWD) 
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GYROSCOPIC CALCULATIONS 
(b) MIDSPAN RESPONSE 
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SEGMENT IN-INIT. IN-FINOL T-flNM. 

(R PM) j <BPM) (BAD) 

1 3000.0 3000.0 6.3 
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FIG. 3.14(a) UNBALANCE RESPONSE WITH INCREASED TRANSLATIONAL 
DAMPING 

(b) MIDSPAN RESPONSE (H-0.05) 
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Fig. 3.16 shows cycles four and five and no indication that numerical 
instability is present in this solution. A small amount of angular 
damping was next considered to control the gyroscopic harmonics and 
prevent them from producing inaccurate results. Five cycles at the 
larger step of 0.05 is shown in Fig. 3.17 and indicates that the 
instabilities have been suppressed. 

For the same Initial conditions as the previous orbit the response 
for ten cycles with an initial step in speed to 10,000 RPNI is shown in 
Fig. 3.18. Fig. 3.19 Indicates the effect of an acceleration rate of 
50,000 RPM/sec. on this same motion for ten cycles. 

These results indicate that the proposed transient solution can 
produce whirl orbits that may be used to study the performance of rotating 
machines. The predicted unstable whirl from the stability analysis 
program was verified with good agreements. The inclusion of gyroscopics 
in a multi-mass system must be handled with extreme care so as to avoid 
numerical instabilities in the solutions and it is suggested that a 
point mass model simulation precede any run including gyroscopics. 

Additional results for flexible supports using this computer code 
are presented in Chapter V. 
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FIG. 3.16(a) UNBALANCE RESPONSE OF OVERHANG FROM THIRD 
CYCLE OF CASE NO. 1215.7101 

(b) MIDSECTION RESPONSE (H=0.01, CYCLES 4-5) 





SEGMENT N-IN1T. N-FINBL T-FINRl 

<RPH) <RPH) <RRP> 

l 3000.0 3000.0 25. 1 
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FIG. 3.17(a) OVERHUNG RESPONSE (GYROSCOPIC DAMPING “ 1IN. 
LB. -SEC. ) 

(b) MIDSECTION RESPONSE (H=0.05, CYCLES 4-8J 
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FIG. 3.19(a) TRANSIENT RESPONSE OF OVERHANG WITH 
ACCELERATION FOR TEN CYCLES 

(b) MIDSECTION RESPONSE (H=0.05, CYCLES 4-13) 

















































CHAPTER IV 


STABILITY OF SOFT MOUNTED JOURNAL BEARINGS 

The favorable increase of stable operating speed range as a result 
of support flexibility and damping has been reported in the I iterature 
(16,74,75,76). A recent analysis by Choudhury (43) has examined the 
effect of simple soft mounts on the stability of the plain cylindrical 
short journal bearing. The present analysis will present the equations 
necessary for the time-transient solution. The resulting computer code 
results are then compared with the stability boundaries presented by 
Choudhury. The influence of unbalance on the horizontal and vertical 
bearing system is also discussed and several transient orbits are present- 
ed in connection with that discussion. 

4.1 Description of the System 

An idealized rigid symmetric rotor supported by identical bearings 
on elastic supports is shown in Fig. 4.1. If the rigid rotor is excited 
only in the cylindrical mode then the equivalent system shown in Fig. 

4.2 may be used to examine the resulting response. A fluid-film bearing 
is shown in Fig. 4.2 with a region of cavitation which must be included 
in the solution to simulate the behavior of actual fluid-film bearings. 

The journal motion may be represented by the coordinates as shown 
in Fig. 4.3. The absolute support motion and the relative journal 
motion will be used to formulate the equations of motion. Unbalance 
response may be expressed as a small eccentricity, e^, of the mass center 
from the geometric center of the journal. The journal orbits will be 


TOO 




FIG. 4.1 SYMMETRIC RIGID ROTOR SUPPORTED BY SIMILAR 

FLUID- FILM JOURNAL BEARINGS ON ELASTIC SUPPORTS 



Y 



FIG. 4.2 CROSS-SECTION OF EQUIVALENT RIGID ROTOR SYSTEM 
SHOWING REGION OF FLUID CAVITATION 
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FIG. 4.3 CROSS-SECTION INDICATING NOMENCLATURE FOR JOURNAL 
DISPLACEMENT SHOWING UNBALANCE OF ROTOR MASS 
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presented as the trajectories of the journal center inside of the 
clearance circle as has been the practice in journal bearing transient 
studies (18,41,60). 


4.2 Equations of Motion for the Transient Solution 

The equations of motion of a rigid rotor supported by journal 
bearings can be written if the expressions for the fluid-film bearing 
forces are known or can be calculated. One approach has been to assume 
that the short-bearing theory is applicable when L/D < I. Booker (77) 
has reported that the short bearing equations compare favorably with 
finite difference solutions. The expressions for the forces can be 
written as (18) 



m a 2 tt (w, + io.)(x sin0 - y cosQ) - 2(x cos0 + y sin 0) cos0 

ifd / -5 J 1 { ..we 

(c - x cos0 - y sine) 3 


0=0 


C4.II] 


where the integral is to be calculated numerically and negative pressures 
set to ambient to model cavitation. 

The coordinate system is shown in Fig. 4.3 where the absolute bearing 
and relative journal displacements have been chosen to describe the 
system behavior. The position vectors are given as: 


P b/o = xjT + yi I [4.2] 

p/°=$»/°+p/n [4.3] 

where 

P j/b = x.T t y j [4.4] 

J J 
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The position vector of the mass center of the journal is given as 


^ g/o = ( Xl + x. + e cos(to.t)) T 
J u J 

+ (yi + Y: + e „ sin (&>.+)) j [4.5] 

Jr* J 

The acceleration of the journal mass center is found by direct 
differentiation of the above expression, 

yJ /o = (xj + x. - e to. sin(io.t)) i 
J V J J 

+ Cy i + y. + e to. cos(to.t)) j [4.6] 

J r* J J 





Ul + Xj - » v i } 

<yi + yj + Vj 


sin(to.t) - e to. cos(to.t)) i 
J V J J 


2 

cos (to. t) - e to. sin(w.t)) j 
J V J J 


[4.7] 


Once the bearing force expression and acceleration components are 
known the equations of motion are easily written by applying Newton's 
Second Law: 




1,2 


[4.8] 


where 

i = I refers to the bearing motion 
i = 2 refers to the journal absolute motion 


Hence for the journal equation: 

m 2 (x 2 t + y 2 j> = l F. t + l F. j [4.9] 

J x J y 
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x-component: 


rMxi + x. - e u>. sin(io.t) - e w. cos(oj.t) 
J P J J V J J 


= F q cos(nWjt) + (Film Force) 


x-d ir. 


y-component : 


*♦ *• 


m 2 (yi + y. + e w. cosCu.t) - e u, sin(w.t) ) 
J P J J P J J 


= F q sin(nWjt) + (Film Force) 


y-dir. 


The bearing equations of motion are as follows: 
x-component : 


• • 

m l x l 55 I Fjj = - k x ><i _ c x xi - (Film Force) 


x-d i r. 


y-component : 


m lYl = I F b = “k Vi - c h - (Film Force) 

y y y 


W x 


| y-dir. 

The expressions for the film forces can be expressed as 
linear stiffness and damping at any given instant. That is. 


F = (Film Force) = -[k x, + k y.+c x.+c y.] 
x x xx j xy'j xx j xy y j 


Fy = (Film Force) 


= -Ck y.+k x.+c y.+c x.] 
y yy y j yx j yy y j yx j J 


where 


k xx = 


9F 

X 

A* 

9F 


xy ay, 


[4.10a] 

[4.10b] 

[4. 1 la] 

[4. 1 lb] 
equ i va I ent 

[4.12a] 

[4.12b] 

[4.13] 

[4.14] 
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yx 


9F 

-JL 

9x . 


yy 


!!l 

8y j ' 


and 


xx 


9F x 

9x7 

J 


xy 


9F 
x 

9 *j 


yx 


yy 


9x7 

J 

A 


9y 


j 


[4.15] 

[4.16] 

[4.17] 

[4.18] 

[4.19] 

[4.20] 


The limits of integration and the force components are determined 
by the steady state equilibrium position (42,43). The equations of 
motion may be made dimensionless by defining the following variables: 


• • 



where c = radial clearance of the journal bearing. 

Rewritting the equations of motion and solving for the acceleration 
terms results in 


A 
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[4.22] 


m 2 k c , F Wi 

Yj — — [■ — — w y ^ - — x — y i — — — y 2 - — 2 ] 
1 mi 1 m2<o, 1 m 2 cw. m 2 cco . z 

J J J J 


F o F X 

X. = -Xx + E cos(T) + 2__ cos(nT) + — 0 

j 1 . p m 2 cw. 4 m 2 ca).2 

J J 


[4.23] 


Y. = -Yi + E sin(T) + — — sin(nT) + 

j p m 2 cu.2 


W 5 


m 2 cu) . z m 2 cto . 2 

t J 


[4.24] 


These equations may be integrated forward in time by any one of 
several methods of integrating first-order differential equations (42). 

This means that to solve the above equations, eight first order equations 
must be integrated. The fluid-film forces can be solved at each time 
step and the solution tracked for a given rotor system. This method of 
solution is essential for determining the dynamic loading under different 
transient motion caused by shock or unbalance forcing functions. This 
approach is not the appropriate method to determine the threshold of 
instability. The characteristic equation can be examined for instability 
by solving for the roots and noting when a positive real part of the 
root, X, exists. Moreover, the method of Routh does not require the 
roots to be solved for but only that an array derived from the coefficients 
of the characteristic equation be formed and examined. The results of 
this type analysis (43) produces stability plots that will be discussed in 
the following section. 

4.3 Horizontal Journal Bearing Stability 

The plain cylindrical bearing has received much attention in the 
literature (78-85). Recent investigations (18,41,69) have presented 
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stability maps that have agreed that the stability boundary Is 

approximated at “/“g - 2.5 for short journal bearings loaded up to an 

eccentricity of 0.80. At higher eccentricities the bearing is stable 

at all speeds. The whirl orbit in Fig. 4.4 represents the motion of 

the journal in the clearance circle for five cycles of motion. The 

journal was released with all initial conditions zero and operating at 

an angular speed of 6500 RPM which is very near but below the stability 

threshold (u>/u> = WS = 2.45). The journal is spiraling in toward the 
9 

equilibrium eccentricity e Q = ES = 0.21 1 . 

If a step increase of speed to 10,500 RPM occurs at the end 

condition of Fig, 4.4 the resulting motion is shown in Fig. 4.5(a). The 

journal whirls out In an orbit that has an average whirl rate of close 

to one-half journal speed as indicated by the small timing marks on the 

orbit path. The maximum force has doubled to 145.6 lb. and occurs at 

8.96 cycles as indicated by the asterisk on the orbit. The motion is 

continued for another five cycles in Fig, 4.5(b) showing the orbit forming 

a limit cycle and whirling at slightly less than 1/2 running speed. The 

"t*h 

maximum force occurs at the end of the l5 Tn cycle of motion and has 
increased to almost four times the static loading. 

The equations in Section 4.2 for the soft mounted journal bearing 
have been examined by Choudhury (43) for stability. If solutions of 
the form X. = Aj e* + are assumed the equations of motion produces a 
determinant given by: 
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FIG . 


HORIZONTAL BALANCED ROTOR 


W. 13301 


N « 6500 RPM 
R * 1. 00 IN. 

L * 1.00 IN. 

C * 5.00 MILS 
TRSMfiX * 1.29 
S » 1.733 

S5 * 0.433 


WT * 1.00 
H * 50 18. 

MU»5 * l.COO RETNS 
FHRX « 64.4 LB. RND 

OCCURS AT 0.53 CTCLE 
W5 * 2.45 
ES « 0.211 



.4 JOURNAL ORBIT OF A BALANCED HORIZONTAL ROTOR ON 
RIGID SUPPORTS (N=6500, W=50, 0=0.005, L/D=%) 
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HORIZONTAL BALANCED ROTOR HORIZONTAL BALANCED ROTOR 

M. I MB! M. 1)9*1 

- 10500 RPN HT - 1.00 N - 10500 RPM NT • 1.00 
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The bearing characteristics are calculated at a given equilibrium 
eccentricity and the resulting characteristic polynomial which is eighth 
order in A was examined by Choudhury using the Routh criterion (86). 

Many stability maps were produced by the analysis and one series Is 
given in Fig. 4. 6-4. 8. The series of curves is for a support to journal 
mass ratio of 0.1 and stiffness ratio K g of 0.01, 0.1, and 1.0. The 
effect of support damping is shown by the family of curves for the 
damping ratio Cg. The length to diameter ratio, L/D, was chosen as 0.5 
for this series of plots. The rigid support stability boundary is shown 
by the dash-dot curve and is as described earlier in this discussion. 

The constant load lines given by the dashed curves represents the 
eccentricity that a particular journal would have as the journal speed 
is varied. 
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FIG. 4.6 STABILITY MAP FOR JOURNAL ON ELASTIC DAMPED 
SUPPORTS (L/D=%, K b =0.01, Mi/M 2 =O.I) 
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ECCENTRICITY, 6o CDIM/D 


FIG. 4.7 STABILITY MAP FOR JOURNAL ON ELASTIC DAMPED 
SUPPORTS (L/D=V. Xb=0.1, Mj /M 2 =0.1) 
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ECCENTRICITY, Co CDIM. 1 ] 


FIG. 4.8 STABILITY MAP FOR JOURNAL ON ELASTIC DAMPED 
SUPPORTS (L/D=>i, K b *1.0, Mi/M 2 =0.1) 
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Examination of this series of plots reveals that there is an opTimum 
support stiffness and damping that will produce the greatest stable 
operating range. The lightest stiffness ratio (Kg = 0.01) appears to 
produce a great increase of oj s even if the support damping is very light. 
Consider the effect of a support stiffness of Kg = 1000.0 I b/ In, damping 
C D = 20.0 lb-sec/in., and weight of 5 lb. Considering these values for 

D 

the example of Fig. 4.4: 

a) = /g/c = 278 rad/sec 

g a 

mi/m2 = 0.1 


r _ 1000 

B " ^°/ 386 a 
c 


0.1 



20 


50 

386 


(278) 


0.556 


The stability boundary should occur at co s = 12 as indicated by the 
stabil ity maps of Fig. 4.6 for Kg = 0.1. For the case of journal speed, 

N = 10,500 RPM, WS = w s = 3.96 and should be stable for this support 
system. Fig. 4.9(a) shows four cycles of the time transient of the 
journal on the support system and it is apparent that the journal has 
been stabilized. The support motion is given in Fig. 4.9(b) and has no 
instability apparent in the motion. The high bearing force is completely 
damped out by the support system. The force to the support is only 6.1 lb. 
over the static loading of the system. 

For a rotor speed of 31,500 RPM the stability parameter o> s = 11.86 
and is at the threshold. The orbit shown in Fig. 4.10(a) indicates that 
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1?7 


FIG. 4.9(a) JOURNAL RELATIVE MOTION SHOWING STABILIZED 
RESPONSE (N=l 0,000, WS=3.96) 

(b) SUPPORT SYSTEM RESPONSE 
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FIG. 4.10(a) ROTOR RELATIVE MOTION AT THRESHOLD OF 
STABILITY (N=31 ,500, WS=11.9) 

(b) SUPPORT RESPONSE AT THRESHOLD FOR FIVE CYCLES 




the journal is indeed unstable and growing at a very slow rate. The 
support motion is negl ibible and is absorbing the large journal forces 
very well but is no doubt increasing at the end of the fourth cycle. 

For a journal speed of 45,000 RPM the system is we 1 1 above the 
threshold and as indicated in Fig. 4. 1 1 (a) the rate of growth is greatly 
increased. The bearing force is 20 times the static loading with the 
support force at the fifth cycle given by 21.9 lb, over the static 
loading. The support motion given in Fig, 4.11(b) is, even in this case, 
still negligible after five cycles of motion. 

Stability analyses such as that conducted by Choudhury must treat 
balanced rotor systems and do not give any indication of the behavior 
under external loading or the degree of stability. Some indication of 
the rate of growth can be obtained If the roots of the characteristic 
equation are solved for In addition to applying the Routh criterion. 

The time transient technique is capable of giving this type information 
and is the reason this type solution is so important to practical design 
theory. Unbalance or any other timewise descri bable force is easily 
included in the solution. Considerable attention has been given to the 
rigid mount short journal bearing time transient behavior (18,41). The 
following discussion will extend this rigid mount journal to include the 
effect of support flexibility on the unbalance response. 

Fig. 4.12 represents five cycles of motion of an unbalanced rotor 

running at 6500 RPM. This Is the same rotor that was examined earlier 

in this section. The journal mass center is eccentric by 0.20 of the 

clearance (ie. EMU = E = 0.20).' The orbit indicates that the motion is 

y 

a combination of synchronous and fractional frequency whirl. The timing 
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HORIZONTAL BALANCED ROTOR « «nt HORIZONTAL BRLANCEO ROTOR *>.w» 

- HSOOO RPH. TRD-B • N - 15000 RPR. TRO-8 

J - 50.0 LB. TRD-S - MJ - 50.0 LB. THO-S • 

a- 5.0 LB. HUoS » 1.0000 RETNS MB- 5.0 LB. HU*5 - 1.0000 RETHS 




OOO’t 499*0 KC'O 000*0- CTC‘0- 499*0- OOO'H 

■‘MO -X 


2 2 
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FIG. 4.11(a) ROTOR RELATIVE MOTION FOR FIVE CYCLES 

SHOWING JOURNAL INSTABILITY (N=45,000, WS=17) 

(b) SUPPORT TRANSIENT FOR FIVE CYCLES (N=45,000) 




HORIZONTAL UNBALANCED ROTOR 


*>. 2J79H 


N * 

6500 RPM 

NT « 1.00 

B « 

1.00 IN. 

H * 50 L8. 

L « 

I. 00 IN. 

ttJmS * 1.000 fONS 

C .• 

5.00 NILS 

FKftX « 135.2 L8. ftNO 

TRSHRX « 2*70 

OCCURS RT 0.86 CYCLE 

S * 

1,733 

NS « 2.45 

ss • 

0.433 

ES * 0.211 

ENU« 0.20 

FU « 59.95 LB. 

su « 

1*446 

FURRTIG « 1.20 

TROMfiX = 2.26 

ESU « 0.244 



marks indicate that the average whirl is somewhat less than synchronous 
and the inner loops indicate the presence of a fractional whirl component. 
As the rotor speed is Increased to 10,500 RPM Fig. 4.13 indicates that 
the fractional whirl dominates the motion as the inner loops disappear 
and the timing marks begin to fall two per whirl cycle. Consider the 
following support system for this example. 

Let 

mass of bearing - 5.0 lb.; mi/m 2 - 0.1 
stiffness of support = 10,000 Ib/in; Kg = 1.0 
damping of support - 200 lb-see/ in; Cg = 5.56 

From the stability map of Fig. 4.8 the system should be stable at 
6500 RPM (co s = 2.45) and possibly close the threshold at 10,500 RPM. 
Figure 4.l4(a)shows the elliptical steady state orbit for the case of 
6500 RPM. The next figure, 4.14(b) indicates the corresponding motion 
of the support system (ie., bearing housing mass). 

When orbits are calculated by any type integration scheme there is 
always the possibility that the results are meaningless. To illustrate 
this fact, the case just considered was run with a larger time step and 
produced the orbit of Fig. 4.14(c). It Is quite obvious that the orbit 
is not correct due to the oscillations starting near the end of the first 
cycle of motion. The next figure, 4.14(d) is not nearly as obvious that 
the solution is incorrect. The support motion is very nearly the same 
as the orbit of Fig. 4.14(b). Extreme caution must be exercised by the 
researcher when producing the time-transient solution for just these 
reasons. 
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HOBIZONTRL UNBRLRNCED ROTOR 


21731 


N « 10500 BPH 

n « 1.00 in. 

L « 1.00 IN. 

C « 5.00 MILS 
TRSMRX » 4.97 
S - 2.800 

5$ » 0.700 

EMU - 0.20 
SU » 0.895 

TROMBX « 1.59 


WT - 1.00 
u m *;n ) n 

MU«5 - 1.000 RETN5 
FMRX » 248.7 LB. RND 

OCCURS RT 8.95 CYCLE 
WS « 3.96 
ES * 0.139 
FU * 156.45 LB. 

FURfiTlO - 3.13 

ESU « 0.342 



FIG. 4.13 JOURNAL ORBIT OF AN UNBALANCED ROTOR ABOVE 

THE STABILITY THRESHOLD FOR CYCLES 6-10(N~10, 500) 
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SUPPORT TRANSIENT MOTION (N=6, 500) 
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FIG. 4.14(c) JOURNAL RELATIVE MOTION SHOWING NUMERICAL 
INSTABILITY FOR SOLUTION TIME STEP TOO LARGE 
(d) SUPPORT RESPONSE FOR NUMERICALLY UNSTABLE 
SOLUTION 




For the support system being considered and journal speed of 10,500 
RPM, Fig. 4.15(a) indicates that the whirl is dominated by a fractional 
whirl and is indeed beyond the threshold speed. The corresponding support 
motion is shown in Fig. 4.15(b). 

Figure 4.16 was an orbit run as a check case while the computer code 
was being developed. The support system was frozen (K = I million, 

C g = 100) and the resulting orbit can be compared to the orbit of Fig. 

4.12 which was presented in (18). 

A lighter support system has been indicated (19) to give better 
steady-state response attenuation. Reducing the stiffness to K g = 1000 
Ib/in and damping to 20 Ib-sec/in, the unbalance response is given by 
Fig. 4.17 a for an unbalance eccentricity of = 0.20. The improved 
response attenuation is readily apparent and the forces transmitted are 
moderate since the unbalance load is 156.6 lb. for this example. The 
support motion and forces are referenced from static conditions and 
Fig. 4. 1 7 (b) i ndicates the bearing housing motion for this particular 
case. The support experiences a maximum force of 33.5 lb. over the 
static loading (55 lb.). 

An increase of the unbalance eccentricity to 0.80 produces an 
unbalance load of 626.4 lb. Fig. 4.l8(a)gives the transient and 
resulting synchronous orbit and a maximum force to the bearing of only 
474.7 lb., less than the unbalance load. The corresponding bearing 
motion is given in Fig. 4. I8(b)and the maximum force transmitted is 
only 136.7 lb. above the static loading. 

This type unbalance attenuation of both amplitude and transmitted 
forces is very desirable when considering the extended life of a machine. 
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FIG. 4.15(a) JOURNAL TRANSIENT RESPONSE FOR OVERDAMPED 
SUPPORT (N=l 0,500, EMU=0.2) 

(b) SUPPORT RESPONSE (KBX=KBP=1 0,000, CBX*= 
CBY=200) 




HORIZONTAL UNBALANCED ROTOR 


H « 6500 RPM. 

UJ • 50.0 LB. 

HB » 50.0 LB. 

EMU - 0.20 
FO « 0.00 LB. 

FHX « 0.00 LB. 

FHT « 0.00 LB. 

KBX.-3 * 1000 LB/ IN 
CBX « 100.0 LB-SEC/ IN 


M0.1U732 

MUs5 « 1.0000 RETNS 

L/D « 0.0 

CL « 5.00 MILS 

EN ■ 0.00 

ENX « 0.00 

ENT » 0.00 

KBTb-3 « 1000 LB/ IN 

CBT - 100.0 LB-SEC/ IN 



FIG. 4.16 JOURNAL TRANSIENT UNBALANCE RESPONSE FOR FROZEN 
SUPPORT SYSTEM (N=6,500„ EMU=0.2) 


128 




HORIZONTAL UNBALANCED ROTOR HORIZONTAL UNBALANCED ROTOR . 0.0.1m 

N = IOSCO RFM. THO-B - 1.531 N - J0500 RPM. TRP-fc ■> 1 -SOI 

KJ » $0.0 LB. TRO-S • 0.21H HJ - 50.0 LO. TfiD-5 - 0.214 

KB « 5.0 LB. MUa5 - 1,0000 RETNS KB - 5.0 LB. MUi.5 • I.OCOO RET NS 

EMU - 6.20 L/0 » 0.500 EMU * 0.20 L/0 - 0.500 



“dd 

z o >- >- 

S£85j 


S' 


o o o o 

^°.do ft 
in o . 

cn 

a u > u 

* 1 >c >- >5 >- { 

dsgsss: 


"V. CO CD 
ZO-I-J 

*— * ll I 

^ to in 

« * • 00 I • « 

co . to cq _j to an <n 
JID J J -J 3* <0 
■ _l CM 

-d°° . 8 , , 

it n ? h co co 

£££E£8££ 



129 


FIG. 4.17(a) JOURNAL UNBALANCE RESPONSE FOR IMPROVED SUPPORT 
SYSTEM (N=l 0,500, EMU=0.2) 

(b) SUPPORT SYSTEM TRANASIENT RESPONSE FOR FIVE 
CYCLES (KBX=KBY=1 ,000, CBX=CBY=20) 




HORIZONTAL UNBALANCED ROTOR HORIZONTAL UNBALANCED ROTOR n.« 

N » 10500 BPM. TRD-B - 0.758 N - 10500 RPH. TRD-0 - 0.7S8 

WJ - 50.0 L8. THO-S - 0.218 !'J - 50.0 LB. TRO-S - 0.210 
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FIG. 4.18(a) JOURNAL UNBALANCE RESPONSE FOR FOUR CYCLES 
(N-l 0,500, EMU =0.8) 

(b) SUPPORT SYSTEM RESPONSE ( KBX=KBY=1 , 000 , 
CBX“CBY=20) 




In many cases the forces transmitted Increase when the amplitudes are 
reduced butthis case demonstrates that damping can be designed into 
nonlinear fluid film bearings to reduce the transmissi bi I ity to the 
bearing surface as well as to the support housing. 

4.4 Vertical Bearing Whirl 

The balanced vertical journal bearing is unstable at all speeds. 

Yet there are machines running in the vertical mode that exhibit 
reasonable response. It has been shown (18) that net external unidirec- 
tional loading can improve the stability of rigidly mounted journal 
bearings. The addition of a damped support system does not improve the 
stability of the vertical balanced journal bearing. The system that 
was stable in the horizontal mode is shown to be unstable in Fig. 4.20(a). 
The journal was given initial conditions 
X = Y = 0.0 
X = -Y = 0. 15 

where 

X = x/c, Y = y/c 
X = x/(cu>), Y = y/(cui) 

The transient response indicates that these initial conditions 
were compatible with the natural response. Although the system is 
unstable the rate of growth is considerably less than that of the rigid 
case. Fig. 4.19. Although the support does not stabilize the journal, 
it does effect the degree of instability. The corresponding support 
response is given in Fig. 4.20(b). It was shown in (18) that an unbalance 
eccentricity greater than 0.16 was required to cancel the nonsynchronous 
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VERTICAL UNBALANCED ROTOR 


Nft, SH 

N • 10500 RPH MT * 0.00 

R • 1.00 IN. H - SOLS. 

L • 1.00 IN. MU»5 - 1.000 BETHS 

C - 5.00 NILS FMftX • 1H8.8 LB. PHD 

TRSHRX « 2.98 OCCURS RT 9. HO CTCLE 

S • 2.800 HS > 3.96 

SS - 0.700 ES • 0.139 

EHU • 0.01 FU • 7.82 LB. 

SU « 17.897 FUBBTIO ■ 0.16 

TRCHRX » 19.02 ESU - 0.023 



FIG. 4.19 JOURNAL ORBIT OF A SLIGHTLY UNBALANCED VERTICAL 
ROTOR FOR 10 CYCLES (N=l 0,500, EMU=0.01) 
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instability whirl and produce a synchronous limit cycle. It was also 
demonstrated that even the unbalance did not alter the fractional whirl 
above the correspond i ng horizontal stability threshold. The following 
series of curves will have an unbalance eccentricity of 0.2. A support 
stiffness of 10,000 Ib/in and damping of 20 I b-sec/in. resu Its in the 
transient orbit shown in Fig. 4.21 (a) for a speed of 6500 RPM. The 
support motion is double that of the journal relative motion as indicated 
in Fig. 4.21(b). If the journal were on rigid supports the motion would 
be stable synchronous as indicated in Fig. 4.22 with the maximum force 
less than the unbalance load. This example illustrates that an improperly 
designed damped support can give undesirable response. If the damping is 
increased to 200 Ib-sec/ln the journal relative motion is shown in Fig. 
4.23(a)and indicates that the amplitude has been reduced by the increased 
support damping. The support motion is also suppressed as indicated in 
Fig. 4.23(b) . 

The preceding analysis of horizontal journal bearings indicated 
that a lighter support stiffness would give improved response. Figure 
4. 24 (a) shows the transient of the improved support system operating at 
10,500 RPM. This lighter support stiffness of 1000 Ib/in gives a 
highly attenuated relative journal response. The support system motion 
given in Fig. 4.24(b) has an initial transient that is damping out to a 
reasonably small whirl orbit. The forces to the support is noted to be 
about one-fifth of the unbalance load. The bearing force is also less 
than the unbalance loading. The rigid support journal is shown in Fig. 
4.25 with a substantially larger response and also a larger maximum 
force transmitted to the bearing support than that of Fig. 4.24(a). 





VERTICAL UNBALANCED ROTOR 
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C * 5.00 NILS 
TRSMRX > 1.09 

S - 1.733 

SS « 0.433 

EMU - 0.20 
5U - 1.446 

TBOIW - 0.31 


•a. visas 

NT - 0.00 
N » 50 LB. 

MU.S « 1.000 PE7NS 

FM3X - 54.6 LB. fWO 

OCCURS RT S.76 CTCLE 
NS - 2.45 
ES « 0.211 
FU « S3. 95 LB. 

FUWTIO - 1.20 

ESU - 0.244 



FIG. 4.22 JOURNAL ORBIT OF AN UNBALANCED VERTICAL ROTOR 
FOR CYCLES 6-10 SHOWING MOTION APPROACHING 
SYNCHRONOUS (N-6,-500, EMU=0.2) 
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FIG. 4.24(a) VERTICAL JOURNAL UNBALANCE RESPONSE FOR 
IMPROVED SUPPORT SYSTEM (N*=l 0,500, EMU=0.2) 
(b) SUPPORT SYSTEM RESPONSE (KBX«KBY=1 , 000, 
CBX«CBY=20) 




VERTICAL UNBALANCED ROTOR 


«L *11583 


H » 10500 RPM 

a « 1.00 in. 

L » 1.00 IN. 

C * 5.00 MILS 
TRSKRX - *1.02 
5 * 2.800 

SS • 0.700 

EHU - 0.20 
SU • 0.895 

TROHBX « 1.28 


«T » 0.00 

H - 50 LB. 

KU*5 * 1.000 REYNS 
FHRX * 201.0 18. mO 

OCCURS RT 8. *19 CYCLE 
MS * 3.96 
ES * 0.139 
FU » 158. *15 LB. 

RJRRTIO » 3.13 

ESU « 0.3H2 



FIG* 4*25 JOURNAL ORBIT OF AN UNBALANCED VERTICAL ROTOR 
FOR CYCLES 5-10 (N«l 0,500, EMU-0.2) 
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An increase of unbalance eccentricity to 0.80 produces an unbalance 
loading of 626 lb. The resulting response is given in Fig. 4.26(a)and 
has a maximum force of 538 lb., less than the unbalance loading." The 
support transient is shown in Fig. 4.26(b)where the maximum force to 
the support housing is only 139 lb. 

The influence of unbalance on vertical bearings has recently been 
found to be effective in producing stable operation of vertical water 
pumps for use in nuclear power plants. This design procedure is just 
the opposite of what is desirable for horizontal machines where unbalance 
induces a synchronous whirling motion. This particular case of vertical 
bearings is another example of the power of the time-transient solution 
in real bearing design. 
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CHAPTER V 


MULT I -MASS FLEXIBLE ROTOR WITH LINEAR 
SOFT MOUNTED BEARINGS 


5.1 Description of the Rotor Model 

A rotor shaft supported by two similar bearings can be reduced to 
the study of motion In a plane If the excitation to the system Is 
symmetrical about the rotor midsection. The continuous rotor mass may 
be lumped at the shaft center and journal locations as illustrated In 
Fig. 5.1. This model then allows the effects of cross coupling Q, 


internal damping Cj, unbalance, and bearing housing flexibility to be 
studied from a relatively simple mathematical model. The coordinates 
used to study the response are given in Fig, 5.2. The bearing housing 


motion is noted as (xj, y i ) , the journal relative motion as (x., y .), 

J J 

and the rotor absolute motion as (xj, y 2 > • 

The influence of the support flexlbi I ity and damping with the 
journal mass neglected was discussed in Section 1.2 and gives a good 
basis for understanding the results to be obtained from this analysis. 


5.2 Equations of Motion 

The system represented in Fig. 5.1 is best described in terms of 
the shaft stiffness and equivalent mass, and m 2 respectively. If 
each journal has an equivalent mass rn. and each support has mass mi, 
then the equations of motion are easily derived (59) and result in 
the following equations: 





FIG. 5.2 CROSS SECTION OF ROTOR MIDSPAN INDICATING 
DEFLECTION NOMENCLATURE 
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Rotor Equations: 


m 2 x 2 + c s ^2 + c j (>'2 ~ Xj - Xj) + k s <x 2 - X! - Xj) 


+ Qy 2 + a)C.(y 2 - yj - y.) - m 2 e w 2 cos(wt) 
■ J 


m 2 y 2 + c s Y 2 + c.(y 2 - yi - Yj > + k s (y 2 - Yi - Yj) 

- Qx 2 - wc. (x 2 - Xi - x.) = m 2 e 0) 2 sin (art) 

* J ^ 


Journal Equations: 

2m. (Xj + x.) + (c. + 2c )x. - c.(x 2 - xj) 

J J * ** J * 

+ (2k xx + k s )Xj - k g (x 2 - x x ) - Cj 0 i(y 2 - yi - 

+ 2(k y . + c y . ) = 0 
xy y j xy y j 


2m j (yi +y.) + (c. + 2c yy )yj - c.(y 2 - y. ) 

+ (2k yy + k s )yj - k s (y 2 - y.i ) + c I w(x 2 - x x - 


+ 2( Vj 


+ c x.) = 0 
yx j 


Support Equations! 

2mjXi + 2m. (xj + x.) + (2cj x + c.)xi - c.(x 2 - x.) 

J J * 1 J 

+ (2k lx + k s )x x - k g (x 2 - Xj) - c.w(y 2 - yx - y.) = 


t5-.ll] 


C5.23 


C 5 . 3 : 


C5.4D 


0 C5.53 



2mm + 2m j .(y 1 + y^) + (2c ly + c,)yx - c f (y 2 - y-j ) 

+ (2kiy + k s )yi - k s Cy 2 - yj) + C|U»(x 2 - x x - x. ) = 0 £5.6] 

The support equations can be reduced to the following form by 
using Eqs. 5. 4-5. 5. 

">i*i + C 1X X! + k lx xj - c^Xj - k^x. - c^Vj - k xyVj . = 0 C5.73 

m ivi + c i V yi + k iyyi - c yy yj - k yy yj - c yx *j - k yx *j = ° Cs-e] 

5.3 Steady State Solution 
5.3.1 Amplitudes of Motion 

The equations of motion of the symmetric three-mass rotor may be 
solved for the steady-state solution by standard procedures available 
in the literature (87). If solutions of the following form are assumed, 

x 2 = A cos(wt) t B sin(a)t) [5.9] 

y 2 = C cos(u>t) t D sin(wt) C5.I0] 

x. = E cos(a)t) t F sln(wt) C5. 1 1] 

vl 

y. = G cos(wt) t H sin(wt) C5.I2] 

J 

Xx = K cos(wt) t L sin(wt) C5.I3] 

yx = N cos(a)t) + P sln(wt) C5.I4] 


the resulting equations will be a 12 x 12 system of equations which may 



be represented In matrix notation as shown on the following page. 

If the following variables are defined, 

ZZ = m 2 U) 

ZZ2 = m 2 (o 2 

Zl = c /ZZ, Z2 = c./ZZ, Z3 = k /ZZ2, F4 = Q/ZZ2 
s I s 

Z5 = m 2 /m . , Z6 = <c. + 2c )/(2 x Z25, Z22 = Z2/2 
J * XX 

11 = (2k + k 1/(2 x ZZ2), Z32 = Z3/2 

XX s 

Z8 = (c. t 2c )/ (2 x Z2) , Z9 = (2k + k )/(2 x ZZ2) 

j yy yy s 

Z 1 9 = k 1722 , Z20 = c 111 , Z2I = k /ZZ2, Z23 = c 111 
xy ' xy ' ' yx ' yx 

ZIO = m 2 /mi, Zll = k lx /ZZ2, ZI2 = c lx /ZZ, ZI3 = k /ZZ2 

XX 

ZI4 = c xx /ZZ, Zl 5 = k ly /ZZ2, Zl 6 = c ly /ZZ 
ZI7 = k yy /ZZ2, ZI8 = c yy /ZZ 

then the elements of the matrix may be obtained from the listing of the 
computer program SSFROLS in Appendix E as the array ACI,vJH. The maximum 
amplitudes in the x and y coordinate directions are then given by 


|x 2 | = X 2 = /A z t B z 


C5.I53 


|y 2 | = Y 2 = /C2 t D2 


C5.I63 


with similar equations for the journal and support equations. 

5.3.2 Phase Angles for Elliptic Orbits 

If a particular motion is described by the equations: 


\kl 


x = A cos(o)t) t B sin (tot) 
y = C cos(uit) t D sin(wt) 


C5.I7] 

C5.I83 



zi le • • • JL • • • e le z t B He 



148 





then the resulting motion Is given (87) by the major and minor el I Ipse 


semi-axis, a and b, and angle of Incl 
where 



with 

n = A 2 + B 2 + C 2 + D 2 
8 = AD - BC 

and 

0 = £ arctan (25/y) 

where 

y = A 2 + B 2 - C 2 - D 2 
5 = AC + BD 

These equations may be applied to the 
displacements to describe the motion, 
to write the x,y components of motion 


nation, 0, as shown in Fig. 5.3, 

[5.19] 

C5.20] 


C5.2I3 

rotor, journal, and support 
An alternate representation is 
as 


x = X cos (ut - 8 ) C5.22] 

A 

y = Y sin(wt - 8y') = Y cos(wt - 8y) C5.23] 

These amplitudes and phase angles can be obtained from experimental 
results as indicated in Fig. 5.4 if the direction of whirl can be 
determined. The phase angles are dependent on whether the whirl is 
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FIG. 



.3 ELLIPTIC ORBIT SHOWING ELLIPSE ANGLE 0 
SEMI-AXIS a, AND MINOR SEMI-AXIS b 


, MAJOR 
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FIG. 5.4 CONSTRUCTION FOR OBTAINING CIRCULAR PHASE 
ANGLES FROM ELLIPTIC WHIRL ORBITS 
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forward or backward. If whirl direction Is known, then the amplitudes 
and phase angles are easily measured with a compass, straight edge, and 
protractor as shown in Fig. 5.4. The inverse problem of trying to 
determine the unbalance level using an orbit trace from an oscilloscope 
requires that these phase angles be determined for an orbit that is 
usually noncircular. Computer codes have been developed (88) which will 
calculate unbalance correction if given amplitudes and phase angles for 
a coordinate direction in two planes. 


5.3.3 Forces Transmitted 

The forces to the bearings and support system can be found with ease 
once the bearing and support amplitudes are computed. The force to the 
bearing is expressed as: 


^Bearing * k xx x j + c xx X j + k xy^j + C xy\j ) 1 

+ lk yy y J + c yy*J + V + V’* 15 - 241 

If the following variables are defined (corresponding to computer code): 

ZI9 = k /(m 2 w- 2 ) 
xy 

Z20 = c /(m 2 w) 

. aY . 

Z2I = k /(m 2 w 2 ) 
yx 

Z23 = c / (m 2 co ) 
yx 

then the dimensionless force expression is 

F B = F Bearing /m2 “ 2 = CZI3 ' x j + 214 * *j + Zl9 * Yj + 220 ' 



This can also be expressed by the component forces in terms of the 
amplitudes of the formulation of page 147/ 

F = (ZJ3 • E + ZI4 • F + ZI9 • G + Z20 • H) cos(wt) 

BX 

+ (ZI3 * F - ZI4 * E + Z 1 9 • H - Z20 • G) sin(a)t) [5.263 


F by = (Z 1 7 • G + ZI8 • H + Z2I • E t Z23 * F) cos(wt) 
+ (ZI7 • H - ZI8 * G + Z2I ’ F - Z23 • E) sin(wt) 
or 

f bx - f ba cos(T) + f bb S,n(T) 

F BY = f bc coslT) + f bo Sln(T) 

where 



/F 2 + p 2 
r BA BB 


[5.273 

[5.28] 

[5.293 

[5.303 



and the major and minor semi-axis can be found by equations similar to 
the equations for amplitudes of motion. 

The support forces may be found by defining the following variables. 
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k lx 

c lx ,, 

k ly 

c ly 


= Zll, 

= Zl 2, 

— — = Zl 5, 

— = Zl 6 


m 2 u 2 


rr^u 2 

n^u 



• 





F $ = (Zll • xi + Z 1 2 • Xj ) I + (ZI5 • y a + Zl 6 • y^j 


C5.32] 


F $x = (Zll • K + ZI2 • L) coswt + (Zl I • L - ZI2 • K) sin tat) 


= F sa coswt + Fgg si nut 


C5.33] 


F $ y - (Zl 5 • N + ZI6 • P) cosut + (Zl 5 • P - Zl 6 * N) sin tat) 

= F^ cosut + Fg D si nut C5.34] 

Hence , 


n S 

Y S 


= F 2 
r SA 


= f 2 

r SA 


+ F 2 + F 2 + F 2 * 
r SB r SC r SD ' 


+ F 2 - (F 2 + F 2 ) • 
r SB ^ SC SD *> 


0 S = F SA F SD " f sb f sc 
h = F SA F SC + F SB F SD 


and the solution may be obtained In terms of the elliptic characteristics 
as before. 

The previous equations for forces give expressions which have units 
of length. A transmisSibi I ity may be defined as the force divided by 
the unbalance loading for a rigid rotor. Then 


TRB - a FB * 2 / E „ 


maximum force to bearing 
unbalance load 


C5.35D 


apg = major semi-axis for the force expression . 
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where 



L ? kew i se 


TRS = a ps x 2/E^ £5.363 

For a linear system the dead weight component may be removed from 
the solutions without effecting the resulting transmissi bi I ity solutions. 
It must be remembered that the system weight must be added to the above 
forces to give the true loading values. 


5.4 Transient Analysis 

5.4.1 Equations for Transient Solution 

The equations of motion as presented in Section 5.3.1 are easily 
solved in terms of acceleration. The linear equations may then be 
solved for the time-transient solution by integration using any of 
several standard procedures available for digital computer simulation. 

The program SMFROLS developed for the transient solution is given in 
Appendix E. Since the equations are linear the solution is obtained 
very quickly as compared to an analysis including nonlinear fluid-film 
bearings. 

The solution may be obtained in terms of mils and dimensionless time 
<03 1) by the following transformation. Let 


X = x x 1 000 ; X = x x 1 000/ oi where 



x x x 


1000 or X = 


d 2 (1000x) 
d?2 


i _ d(IOOOx) 

X ‘ dT ’ 


T - ut 


with similar expressions for the y-coordinate variables. Then the 
equations are: 
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*2 - % - 10*) cos(T) - (J^X* - *1 - *j> 


' ( 53 f )(x » ' x > ' x j’ " < ^r >Y * ' ( S ^ ><Y2 * Y * ■ V C5 ' 37] 


Y 2 - Ce y * I0 3 > sin(T) - - <Sp?>^2 " V" *,> 

- 4» y ! - Y > - y j> + ( 4^ >(X2> + & <x * ■ Xi * v C5 - 3811 


.. ^2 C, + 2C . C. . 

X J - - x > - ( SJ>< ( -W^ )X J • W' X2 • Xl> 


2K + K 


<j " ( afer ><x - 2 - Xl> " ( 2i^ KY2 


+ (-2I_)Y. + C^OY.} 
j m2W J 


Y. « -Yi - (— 
j m. 


m 2 G, 4* 2C 

) f (, 1 J 

m. 2m 2 w 


0 - w (Y * - Y >> 


2K + K K C. 

+ Hfer*» Y j - ( ak )(Y2 - Y >> + ( ^ )<X2 

+ ( S^T )X J + } 


m 2 K lX 


Clx. x 


*1 *^57 ){( iiS5- )Xl + ( ii^ >Xl 


<-25*>X. - <^*)X, 


m2W J 


- (-22L)t. - (r^)Y,> 

n)2“ z J m 2 14 J 
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IT12 Kly Cjy 

Yi =-(—){ (-T^-JYi + (r-f-)fi 
mj ni2W^ * m2<o 

- (^%)X. - 0*>X.} 
m2<*) z J rr^w J 


K 


- (JO-w - (J^)Y. 

m^o)* j m2® J 


15.423 


5.4.2 Time Transient Whirl Orbits 

Consider the following example rotor system. 
Example 5.1 

Rotor Weight = W 2 = 675 lb. 

Journal Weight = Wj = 312 lb. 

Support Weight = W \ - 50 lb. 

Bearing Characteristics: 


•S*- 

1,287,000 lb/ in. 


1,428,000 1 b/ in. 

O 

II 

1,200 lb. -sec. /In. 

O 

-< 

-< 

II 

1,290 lb. -sec./ in. 


Support Character I sties: 

K lx ■" Kiy = 130,000 lb./ in. 

C lx - Ciy - 150 lb. -sec. /In. 

The shaft stiffness is 280,000 lb. /in. for this example. The anti- 
symmetrlcal bearing characteristics given are typical fqr a rotor of 
this weight. The critical speeds for rigid support housings are 3630 
RPM and 3647 RPM for the x and y coordinate directions respectively. 

A series of orbits representing steady-state solutions were calculated 
using the transient solution computer code and the steady-state 
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starting conditions as calculated by the program SSFROLS. The mark on 
the orbits represent a timing mark such as that produced by a keyphasor 
probe detecting a mark on the rotor shaft. Examination of the three 
orbits (ie., rotor, journal, and support motion) indicate ; phase angles 
of approximately 30° at the speed of I80(J RPM (Fig. 5.5). The next 
series of orbits at 3600 RPM indicate phase shifts of nearly 180° which 
indicate that the rotor system including the support has passed through 
its first critical speed (see Fig. 5.6(a), (b),(c)). The last series 
of steady-state orbits (Fig. 5.7) at 18000 RPM indicate that the absolute 
rotor motion phase remains at 180°. The journal relative motion phase 
Is approximately 280° while the absolute motion phase is close to 350°. 

The support motion has a phase shift of 360° at this speed of 18,000 
RPM. It is noted that the rotor amplitude is approximately equal to 
the unbalance eccentricity, EU = e^ = 5 mils. 

The maximum force encountered during the corresponding cycles of 
motion has been indicated on the previous series of figures. The 
absolute rotor motion orbit has the unbalance force, FU indicated. The 
force transmitted to the bearing is indicated on the journal motion 
orbits as TRDB. The time of the maximum force is indicated in parentheses 
following the value of TRDB. The maximum force transmitted to the 
support Is likewise indicated on the support motion plots. For reference 
the TRDB values for 1800, 3600, and 18,000 RPM are 2.88, 0.920, and 
0.007 while the TRDS values are 2.984, 1.043, and 0.005 respectively. 

if the support system is rigid then the rotor motion for EU = 0.5 
mils. Cl = 20 Ib-sec/in, and N = 16000 RPM is shown in Fig. 5.8(a) for 
ten cycles of motion and continued for cycles 10 - 20 in Fig. 5.8(b). 
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FIG. 5.5(a) STEADY-STATE ORBIT FOR ROTOR OPERATING AT ONE-HALF THE RIGID 
SUPPORT CRITICAL SPEED 

(b) JOURNAL MOTION FOR FIVE CYCLES 

(c) SUPPORT MOTION INDICATING A PHASE OF 30° 
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FIG. 5.6(a) STEADY-STATE ORBIT FOR ROTOR OPERATING AT APPROXIMATELY 
THE RIGID SUPPORT CRITICAL SPEED 

(b) JOURNAL STEADY-STATE ORBIT SHOWING PHASE ANGLE APPROACHING 180° 

(c) SUPPORT MOTION SHOWING PHASE OF 180° 
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FIG. 5.7(a) STEADY-STATE ORBIT FOR ROTOR OPERATING WELL ABOVE RIGID SUPPORT 
CRITICAL SPEED (PHASER 80°) 

(b) JOURNAL RELATIVE MOTION (PHASE=275°) AND ABSOLUTE MOTION 
(PHASE— 355°) 

(c) SUPPORT MOTION SHOWING PHASE OF 360° 
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FIG, 5.8(a) RIGID SUPPORT ROTOR MOTION SHOWING UNBALANCE 
RESPONSE WITH LIGHT INTERNAL DAMPING (CYCLES 0-10) 
(b) RIGID SUPPORT ROTOR MOTION CONTINUED FOR AN 
ADDITIONAL TEN CYCLES 



The system is highly unstable under these conditions. The journal motion 
given in Fig. 5.9(a) and 5.9(b) indicate the same instability. The 
maximum force for the motion is given as 0.35 of the unbalance loading 
to the bearing. 

If the rotor system in the previous rigid support orbits has a 
damped flexible support of stiffness 130,000 lb/ in and damping of 
300.0 Ib-sec/in the results are shown in Fig. 5.10. The motion indicates 
that the support has absorbed the instability and would reduce to 
circular synchronous motion after the initial transients have decayed. 

Fig. 5.11 represent the system motion for a reduced support damping 
and increased internal damping Cl = 132 Ib-sec/in. The unbalance level 
was reduced and the motion was started from steady-state conditions by 
an initial velocity perturbation on the rotor. The instability is 
observed to begin in the support system and would eventually drive the 
absolute rotor motion to an unstable whirl orbit. 

The computer code SMFROLS can be used as an additional check on 
the program MODELJ developed in Chapter II. The same rotor system as 
Example 5.1 will be used to compare the results of the two programs. 

For a rotor speed of 10,000 RPM and an unbalance of 0.500 mils the 
rotor absolute motion is shown in Fig. 5.12(a) for ten cycles of motion. 
The unbalance load in 958.73 lb. as indicated on the figure heading. 

The corresponding journal motion is given in Fig. 5.12(b) where the maxi- 
mum force transmitted to the bearing support is indicated to be 314 lb. 
This is calculated from the dynamic transmissibi I Ity factor as follows: 

F max = TRDB x FU/2 = 0,656 x 958 - 473 / 2 * 0 = 314 lb. 
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FIG. 5.9(a) JOURNAL MOTION FOR UNBALANCE RESPONSE WITH 
LIGHT DAMPING AND RIGID SUPPORT SYSTEM 
(b) JOURNAL MOTION CONTINUED FOR AN ADDITIONAL 
TEN CYCLES 
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INTERNAL DAMPING AND MOUNTED ON FLEXIBLE, DAMPED SUPPORTS 
JOURNAL RELATIVE AND ABSOLUTE RESPONSE 

SUPPORT TRANSIENT RESPONSE FOR SUDDENLY APPLIED UNBALANCE 
LOAD (N«16,000) 
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FIG. 5.11(a) ROTOR ABSOLUTE MOTION FROM STEADY-STATE RESPONSE EXCITED BY 
A 1.1 IMPULSE LOAD 

(b) JOURNAL ABSOLUTE MOTION INDICATING UNSTABLE SUPPORT SYSTEM 

(c) SUPPORT TRANSIENT RESPONSE TO IMPULSE LOAD ON ROTOR (N=18,000) 
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FIG. 5.12(a) ABSOLUTE ROTOR RESPONSE TO SUDDENLY APPLIED 

UNBALANCE LOAD (RIGID SUPPORT, 5.40Z-IN UNBALANCE) 
(b) JOURNAL TRANSIENT FOR TEN CYCLES (N=l 0,000) 



The force occurs at 0.76 cycles of motion as indicated in parentheses 
fol lowing the TRDB factor. 

The same rotor system was programmed for the multi-mass flexible 
rotor program and the results of the calculated orbits are given in 
Figs. 5.13(a) and 5.13(b). for five cycles of motion. The orbits are 
similar to Figs. 5.12(a) and 5.12(b) for the first five cycles of motion. 
The maximum force to the bearings is predicted by MODELS to occur during 

i 

the first cycle of motion and is given as 290 lb. for each bearing (see 
Table 5.1). This represents an 8 % derivation in the two program on the 
maximum forces calculated. This is due to the fact that the 0.5 mil 
eccentricity is equivalent to 5.4 oz-in unbalance. The orbits of 
Figs. 5.13(a) and 5.13(b) were produced from an unbalance level of only 
5 oz-in and explains the slight difference in the orbits and forces. 

A similar comparison is now made for the flexible, damped support 
system having support stiffness of 130,000 Ib/in and 300.0 Ib-sec/in 
damping. The absolute rotor, journal, and support motion calculated 
by SMFR0LS are given in Figs. 5.14(a), (b), (c). The dashed lines in 
Fig. 5.14(b) represents the journal, absolute motion while the solid line 
is the relative motion. The maximum bearing force to each bearing is 

F B „ = TRDB x FU/2 = 958.473 x 0.402/2 = 193 lb. 
max 

The support force is 

F = TRDS x FU/2 = 958.473 x 0.429/2 = 206 lb. 

max 

The orbits produced by the program MODELJ with unbalance level of 5.4 
oz-in. are given in Fig. 5.15 and are in good agreement with the orbits 
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FIG. 5.13(a) TRANSIENT ROTOR RESPONSE FOR FIVE CYCLES BY 
PROGRAM MODELJ (50Z-IN UNBALANCE) 

(b) JOURNAL TRANSIENT RESPONSE TO IMPACT 
UNBALANCE LOAD OF 5 OZ-IN ON ROTOR 
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TABLE 5.1 ROTOR SPECIFICATIONS FOR UNBALANCE RESPONSE 
ON RIGID SUPPORT BEARINGS INDICATING MAXIMUM 
FORCES TO BEARING HOUSING 





ABSOLUTE ROTOR MOTION JOURNAL MOTION SUPPORT MOTION 
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FIG. 5.14(a) ABSOLUTE ROTOR TRANSIENT RESPONSE TO IMPACT UNBALANCE LOAD 
FOR DAMPED ELASTIC SUPPORT SYSTEM (5.40Z-IN UNBALANCE) 

(b) JOURNAL TRANSIENT ABSOLUTE MOTION FOR TEN CYCLES 

(c) SUPPORT SYSTEM TRANSIENT RESPONSE FOR TEN CYCLES (N=l 0,000) 
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FIG. 5.15(a) ROTOR TRANSIENT RESPONSE WITH 5.40Z-IN UNBALANCE IMPACT LOAD 

(b) JOURNAL RESPONSE FOR FIVE CYCLES 

(c) SUPPORT TRANSIENT MOTION FOR FIYE CYCLES (N=l 0,000) 
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of SMFROLS. Table 5.2 indicates very good agreement with the maximum 
forces calculated and indicates the maximum force to each bearing to be 
190 lb. and the maximum force to the supports is indicated to be 210 lb. 
This excellent agreement of the two programs indicates that the theory 
of Chapter il produces results which are accurate and acceptable for 
study of complex rotor support systems. 

5.5 Stability Analysis of the Three Mass Rotor With Aerodynamic 

Cross Coupling and Internal Friction 

If tf > rotor analyst is primarily concerned with the design of 
flexible support systems to promote stability and for the attenuation 
of the rotor amplitude and forces transmitted due to unbalance, it would 
be difficult to utilize the transient orbit program for the solution of 
the optimum support parameters,' For example, in the case of the rotor 
in Example 5.1 with internal friction it would be desirable to know the 
range of the permissible values of the support system stiffness and 
damping that may be incorporated for optimum performance. It is difficult 
to determine from the observation of the orbits of Figs. 5.10(a) and 
5.11(a) for support damping values of 300 and 150 lb. -sec. /in. respectively 
whether the optimum damping is higher or lower or whether the support 
stiffness should be reduced or increased. It would be extremely expensive 
and time consuming to develope the necessary series of rotor orbits for 
all the ranges of variables required for the desired results. 

For this reason it is therefor desirable to first develop stability 
maps showing the relationship of the various variables on rotor stability. 
Variables of interest would consist of support housing damping and 
flexibility, bearing characteristics, rotor shaft flexibility. Internal 
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damping, and aerodynamic effects. The stability analysis for the linear 
system is identical to the analysis that will be presented in Chapter VI 
for the stability study of the nonlinear squeeze damper. The computer 
code SDSTB of Appendix F may be programmed for given bearing and support 
characteristics and hence may be used to develop© stability charts for 
the given rotor bearing system. 

For example consider the effect of the damped support on the 
suppression of aerodynamic instability as shown in Fig. 5.16 which has 
the value of the real part of the major root plotted versus the support 
damping for various values of support stiffness. The support mass for 
this example is noted to be 15 lb. with the other variables remaining 
the same as given for Example 5.1. This plot is for a rotor speed of 
10,000 RPM and an aerodynamic cross-coup I i ng of 20x000 lb./ in. The 
figure indicates that there is a definite value of support damping that 
will promote stability for a given support stiffness. For example 
a damping of 500 l,b-sec/in. would be optimum for the support stiffness 
of 50,000 I b/ in. As the support stiffness increases fhe required optimum 
damping also increases and the system becomes less stable for the 
Increased support stiffness. If the support system loses the damping 
then the rotor system is easily driven unstable as indicated on the 
left-most portion of the stability map. 

The importance of the tuned support damping is even more evident 
for the increased level of cross-coupling given in Fig. 5.17. It is 
noted that the range of damping that will promote stability has been 
reduced considerably from a range of 100 to 10,000 for the lesser 
cross-coupling value to only a range of 370 to 1600 for the higher 
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STABILITY OF A FLEXIBLE ROTOR WITH AERODYNAMIC 

CROSS COUPLING 
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level of aerodynamic excitation. The stiffer support systems are 
indicated to be unstable regardless of the damping chosen for the 
support system for the case of the higher aerodynamic cross-coupling. 

At this point in the analysis the transient program could be used 
to establish the effect of unbalance and external shock loads to the 
support system. The transient analysis becomes of even greater importance 
when nonlinear system are being analyzed for stability by linear 
approximations. These nonlinear systems may exhibit stability for small 
pertubations but may also give large sustained whirl orbits for suddenly 
applied unbalance and Impact loads. The nonlinear damper will be 
discussed in greater detail in the following Chapter. 
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CHAPTER VI 


MULT I -MASS FLEXIBLE ROTOR WITH NONLINEAR 
SQUEEZE DAMPER SUPPORTS 

This chapter will Introduce the nonlinear squeeze damper support 
system into the model analyzed in Chapter IV with the added feature of 
allowing for rotor shaft flexibility. Thus the model of this chapter 
will be closer to real life machine performance than the linear model 
of the previous chapter. The inclusion of nonlinear elements in the 
transient analysis can be very time consuming unless approximate 
stabil ity boundaries are known. A computer code to determine the 
stability of the rotor model described in Section 6.1 was developed 
and will be discussed in connection with a time transient program for 
the same model. 


6.1 Explanation of the Rotor Bearing Model 

The rotor-bearing system shown in Fig. 6.1 represents a three 
mass symmetric flexible rotor supported by journal bearings suspended 
in squeeze damper housings. The analysis will not consider the 
gyroscopic action of the rotor disk and hence the problem may be 
reduced to the study of motion in a plane. The characteristics of the 
fluid-film bearing have been discussed in Chapter IV and will be 
used to model the bearings in the present analysis. The squeeze 
damper supports will be derived from the same analysis. Fig. 6.2 
indicates the cavitated bearing model and the fluid film damper recess. 
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FIG. 6.1 ROTOR SYSTEM SUPPORTED ON FLUID- FILM JOURNAL BEARINGS 
WITH BEARING HOUSINGS SUSPENDED IN FLUID-DAMPER SUPPORTS 
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FIG, 6.2 CROSS-SECTION OF THE DAMPER STATION SHOWING JOURNAL 

CAVITATION AND RETAINER STIFFNESS IN THE BEARING SUPPORT 
HOUSING 


The bearing ring has an anti-rotation pin and is suspended by assumed 
linear retainer springs which may have unequal stiffness characteristics 
for the two coordinate directions, 

6.2 Equations of Motion 

The equations of motion of a symmetric two bearing rotor may be 
written in terms of the following coordinates. 

x 2 ,y 2 “ absolute displacement of rotor mass, m 2 

Xj,yj - absolute journal displacement 
x b’^b ” bearing bush motion 

The equations are easily derived considering linearized stiffness and 
damping terms to model the hydrodynamic bearings and squeeze bush. A 
simple force balance and inclusion of the internal friction components 
and cross coupling expressions at the rotor station produces the 
f o 1 1 ow i ng equat i ons . 

Rotor Mass: 


m 2 x 2 t k g (x 2 - x.) t lC(x 2 - Xj) t a)lC(y 2 - y j ) 



m 2 e^w 2 cos(wt t $ 2 ) 


( 6 . 1 ) 


•* • 
m 2 y 2 + k g (y 2 “ Yj) + lC(y 2 - yj) - wIC(x 2 - Xj) 


+ c s Y 2 “ Ox 2 - m 2 e^u) 2 sin(tot + <j> 2 ) 


( 6 . 2 ) 
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Journal Bearings: 

• • 

2m J X J ” k s (x 2 “ Xj) - JC(x 2 - Xj) - wlC(y 2 - yj) 

+ 2 V X J ' V + 2k xy (y J ' y b’ + 2c xx (i J ' V 

+ 2c xy <y J ’ V = 0 C6 ' 3] 

2m j y j " k s (y2 " y j J ~ IC(y 2 " yji + <*>ic(x 2 * Xj) 

+ 2k yy ly J " V + 2k yx <x J " x b> + 2c yy<Vj * V 

+ 2= yx <y, - y b » * 0 Cs - 4:l 

Bush Damper 
*« 

2m b x b + 2k lxx x b + 2k lxx x b + 2l Vb + 2c lxx x b + 2c lxy*b 

• 2k xx (x J ' x b> ‘ 2k xy^J - V - 2c xx'*J - V 

- 2c xy <y, - y b > - 0 Cb.53 

2 Vb + 2k lyy y b + 2k lyx x b + 2c lyy y b + 2c ly x x b 

- 2k yy Cy J - y b> - 2k yx U J ' V ' 2c yy ly J ~ V 

- 2c yx (Xj - x b ) =0 [6.6] 

These equations are basically the same as those presented in the 
previous chapter. The difference Is due to the absolute coordinate 
chosen to describe the journal motion. The methods of obtaining the 
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bearing characteristics have been discussed In Chapter IV and the 
squeeze damper characteristics will be presented in the following section. 


6.3 Stabi I ity Analysis 

Manufacturers are presently installing floating squeeze damper 
supports into their high speed machines. Tanaka and Horl (89) have 
recently reported on the stability of the spinning floating bush 
damper. They included rotor shaft flexibility but excluded the journal 
mass from the analysis. Several curves were presented showing that regions 
of instability could be passed through as rotor speed increased. 
Experimental results were presented and the agreement with the analytical 
results was satisfactory. 

The major conclusions of Tanaka-Hori may be summarized as follows: 

(1) Floating bush damper bearings are superior to cylindrical 

bearings due to the damping effect in the bush. 

(2) Stability is generally improved when 

a) the shaft is stiffen 

b) large floating bush radius 

c) large bush clearance 

d) small L/D ratios 

e) small oil supply pressure 1o allow cavitation 

(3) The stability criterion is complicated and a particular 

case must be checked to assure best design. 

As early as 1958 Hill (90) reported on the improved performance 
of small gas turbines running at 36,000 RPM by using slipper bearings 
which actually introduced a damped support system into the bearing 

m 



configuration. Dworskie (91) later reported successful use of floating 
sleave bearings in a 44,000 RPM gas generator. 

In 1968 Orcutt and Ng (92) presented steady-state and dynamic 
properties of plain cylindrical floating ring bearings but the 
equations were reduced by the assumption of synchronous motion which 
limits the validity of these results since fractional frequency whirl 
is to be expected in the fluid-film bearings. Experimental tests 
indicated several different types of instabilities were possible. The 
journal was usually the source of instability but at times the outer 
film was unstable while the inner film was stable. One conclusion 
reached was that the amplitude of the instability was more important 
than a calculated stability threshold speed since low-level whirling 
could be tolerated in any given application. 

The present analysis will consider the damper bearing as a non- 
spinning journal deriving damping only from the fluid film squeeze 
action. Considering a pure squeeze action for a velocity in the +y 
direction the positive film extent is from 0 to t as shown in Fig. 6.3(a). 
For a positive x velocity the positive film extends from-ir/2 to ir/2 as 
indicated in Fig. 6.3(b). The expression for the force is obtained 

from Eq. 4.1 for the case of w. = 0.0. 

J 


F (s) yRL 3 
x,y 2 


J 11 ~2(xcos 8 + ysin6) 

0 (c - xcosQ - ysinS) 3 




C6.7D 


Differentation of this equation produces the damping coefficients 
as fol lows: 
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aF x (x,y = 0) 


'ixx ax" 


_ yRL 3 


tt/2 2 cos 2 0 

J 


—it/ 2 (c - xcosO - ysin0)' 


d 0 [ 6 . 8 ] 


lyx 


3F (x,y = 0) □, 3 tt/2 2cos 0 sin 0 

=a r- I, r« C6.9D 

-ir/2 (c - xcos 0 - ysin 0) 3 


_ _ 8F x^ X _ yRL 3 J 2cos0 sin0 

j 


lxy 9y 


o (c - xcos0 - ysin0) 3 


d0 [6.10] 


c = - 9F y (x ~ °’Y ) = yRL 3 j 2sin 2 0d0 

^ ^ o (c - xcos0 - ysin0) 3 


[ 6 . 11 ] 


It is of interest to note that McGrew (93) has derived a damping 
expression for circular synchronous precession (short bearing solution 
using rotating coordinates) wnich is expressed as 


C = r ft - i 

° ( | _ £ 2 ) 3 / 2 
O 


[ 6 . 12 ] 


Also by assuming e = eto cos (tot) a pseudo stiffness term may be derived 
and is given by 


^ _ 2ywL 3 R 
K o 



[6.13] 


The results of the computer calculated damping coefficients are 
presented in Fig. 6.4(a) along with points calculated from the expression 
derived from McGrew’s analysis. The circular synchronous damping values 
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DAMPING [DIM] 



0.0 0.2 0.4 0.6 0.6 1.0 


€ , ECCENTRICITY [DIM] 


FIG. 6.4(a) FLOATING BUSK DAMPING CHARACTERISTICS VERSUS 
ECCENTRICITY (X=0.0, Y=e) 
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are observed to equate to c lxx while the stiffness term Is actually a 
cross coupled damping force as It must be since the damper can In no 
Instance develop a stiffness In the absolute sense. 

Figures 6.4(b) and 6.4(c) Indicate that preloading the damper 
produces damping expressions that are converging to equivalent values 
for the two coordinate directions. This characteristic may be very 
desirable when trying to design a tuned support system for optimum 
performance. 

Figure 6.5 Illustrates the effect of the additional pseudo stiff- 
ness term on the rigid rotor supported by a damper bearing with re- 
tainer spring rates of 5772 I b/ In. The initial condition was given for 
the zero speed sag of the retainers. At the running speed of 9500 RPM 
the journal lifts up from the steady^state balanced equilibrium 
position and has a stable whirl orbit due to the unbalance level of 
0.08 of the clearance. This lift of the nonlinear fluid damper must 
be considered when designing real machines which do have unbalance. 

The additional psudo stiffness (damping) contribution could shift the 
operating conditions with unbalance into a region where aerodynamic 
excitation could drive the system and produce large unstable whirl 
orbits. 

j 

Once the damping expressions for the floating bush and journal 
characteristics are known the stability analysis Is easily programmed 
on a digital computer. 

Assuming solutions of the form 

x f - A f e X+ ; 1 = 1,6 [6.141] 
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FIG. 6.4(b) FLOATING BUSH DAMPING CHARACTERISTICS VERSUS 
ECCENTRICITY (X=0.10) 
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FIG. 6.4(c) FLOATING BUSH DAMPING CHARACTERISTICS VERSUS 
ECCENTRICITY (X=0.20) 
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SQUEEZE FILM BEARING 

HORIZONTAL UNBALANCED ROTOR 


NO. 1197]} 


N « 9500 RPM 
R = 1.20 IN. 

L = 1.00 IN. 

C - 9.00 MILS 
TRSMAX = 1.14 

KRX = 5772 LB/ IN 
EMU » 0.08 
SU * 1.332 

TRDHAX = 0.62 


WT - 1.00 
M * 28 LB. 

MUa5 « 1.000 REYNS 

FMAX * 31 .4 LB. AND 

OCCURS AT 1.69 CYCLE 
KRY = 5772 LB/ IN 
FU * SO. 71 LB. 

FURATIO = 1.84 

ESU = 0.200 



FIG. 6.5 SQUEEZE DAMPER UNBALANCE RESPONSE FOR RIGID ROTOR 

AND RIGID BEARINGS (N=9,500RPM, W=28L3, C=9MILS, EMU=0.08) 
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where 

i = I corresponds to x 2 

1 = 2 corresponds to y 2 

1 = 3 corresponds to x^ 

i = 4 corresponds to yj 

i =5 corresponds to x^ 

1 = 6 corresponds to y fa 

The homogeneous equations of motion may be examined for stability by 
determining the solutions of the characteristic equation and examining 
the real part for a positive value. The characterist ic equation can be 
obtained by expanding the determinant given on the following page. 

A computer code SOSTB has been developed to carry out the 
necessary operations to determine the stability of a given rotor-bear ing 
system. A listing and explanation of the required input is given in 
Appendix F. 

To demonstrate the capability of the program consider the following 
example case study. 

Example 6. I 


Rotor weight 

- 40 lb. 

Journal weight 

- 5 lb. (each) 

Bush weight 

- 2.5 lb. (each) 

Journal clearance 

~ 5 mi Is 

Journal radius 

- 1 in. 

Journal length 

~ 1 in. 

V i scos i ty 

- 1 x I0~ 5 reyns 

Radius of bush 

- 1.2 in. 


193 



1 

o 

1 

o 

o 


O 

o 

’ 1 
O 

II 

1 

r-t 

CM 

cn 


J* 

cn 

1 

u> 

< 

1 

<C 

< 


< 

< 

< 

1 


>• > 
X X 
3 «X 


> >* 

> > 

Y f 


>- X 

as 

4 ? + >- 

x" 

■& ? 


>» >• 
>* >* 
H IH 

4? f 

>* > 
CM > 

r< O -X 

jQ ^ ^ 

E + + 


r< 

8 8 
V f 


X 

>- 


-X 

I 


"8 8 

ri iH 

* x + x 

% o** x 

Xl w 

E + + 


r< >-> 

~x £ 

•*«> 
& + 


O 


o 


<< 

I 


IT* 

o + 


f< 

CN 

CM N 

>*. «-s 

O O 

M »— I 

4 * 3 

>- I 

CM >- (0 

* O -* J 

E + + + 


*< 

^3? 

V f 


>* >. 
> >. 

Y f 


a « 

*< .x 

« i 


o 

*-i 

3 


r< 

CM 

CM \ 

c5 o 

*H M 

+ x| 

CM X U) X 
«■< O -X X 

T^* w Jd 

E + + + 


r< X 
X >• 

>-.x 

u + 


-X 

I 


X X 
>• >- 
o -X 


o 


CM O 
*< *~« ID 

CM -X 

£ *f + 


CM 

O 


CM 

O 


-X 

I 


CM 


</> 

« n 

r< »-< <fl 

CM w ^ 

l!- + . !„ 


o 

►-« 

3 

dv 


O 

CM j3 

O 

-X 

i i 


CM 

X 


<M 

> 


-O 

X 


x> 

>* 


19 ^ 





k ix = k iy = re ^ a i ner stiffness - 5772 lb/ In. 

The bush clearance was varied from 5 to 30 mils for each of two 
shaft stiffness values of 5,000 and 500,000 1 b/ in. These produce 6/c 
values of 1.6 and 0.016 respectively. The result from Ruhl's analysis 
(40) of the flexible rotor on short journal bearings is shown in Fig. 

6.6 for reference. Comparing the above values of 6/c to the flexible 
rotor stability chart of Ruhl the stability of the flexible rotor (6/c = 
1.6) should be lowered considerably from that of a rigid rotor. 

The results of the analysis is given in Fig. 6.7 where the stability 

boundaries for both shaft stiffness values are plotted. For a bush 

clearance of 5 mils the stability of the flexible shaft is indeed. 

considerably reduced from the rigid rotor values. This is as predicted 

by Ruhl. As the bush clearance is opened up the stability boundary 

increases first for the flexible shaft and then for the rigid shaft. 

For a clearance of 12 mils the stability boundary has increased from 
a> 

/g7c7 = (1.0, 2.54) to (8.0, 13,0) respectively for 6/c = (1.6, 0.016). 

J 

As the clearance increases their is a lower region of instability which 
can be passed through to a region of stable operation. As the rotor 
speed increases the whirl reduces from 0.50 to 0.25 at 17 mils and 
encounters another instability of whirl frequency of approximately 0.50 
as predicted by the stability program. As the clearance is opened 
further the lower instability region increases in width and finally 
converges with the upper instability boundary at higher bush clearance 
values (~ 28 mils for 6/c = 1.6). This plot was developed by keeping 
the retainer spring rate fixed and increasing the bush clearance. This 
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of course produces results which vary from an analysis that has the 
bearing always centered in the bush clearance. 

This phenomenon of passing through regions of instabi I ity with the 
damped supports has been reported in the I iterature for experimental 
results (54,74) and increases the significance of this analytical simula- 
tion results. The results of a time transient computer code to verify 
certain boundaries shown in Fig. 6.7 will be presented and discussed In 
the following section. 

6.4 Transient Behavior of a Rotor Supported in a Floating Bush Damper 

The equations of motion including the nonlinear journal bearings 

and squeeze bush damper force evaluations have been programmed for 

verification of the stability program SDSTB. This computer code, TMASSNL, 

which calculates and plots the transient solution is given in Appendix G. 

The equations programmed were written in terms of journal relative motion 

for the transient solution and are identical to Eqs. 5.37-5.42. The 

~f"h 

program has an option to use either a 4 n order Runge-Kutta or an Euler 
integration technique. The type integration chosen determines the step 
Size necessary for a realistic whirl orbit. The step size required 
depends on the loading, rotor size (weight) and stiffness, and the type 
support chosen for the simulation. As the number of degrees of freedom 
Increases the step size generally must be decreased to insure a valid 
solution. 

Several cases have been computed that were chosen to verify the 
stability map of Fig. 6.7. Consider first the rigid rotor (<S/c = 0.016) 
with a bush clearance of 17.5 mils. From Fig. 6.7 the system upper 
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stability bound is at a> s = 12.3 or 32,000 RPM for our system. Fig. 

6.8(a), (b), (c) indicates two cycles of motion at 29,000 RPM of the 
absolute rotor, journal relative, and support motion respectively. The 
system was started near the steady-state equilibrium position and the 
journal relative motion is approximately a half-frequency whirl and the 
damper is moving in a non-circular path as indicated in Fig. 6.8(c). The 
next series of figures continue the motion for cycles 2-5 and indicates 
that the rotor motion Is negligible (Fig. 6.8(d)) while the journal 
relative motion is still whirling but decreasing as indicated by the orbit 
of Fig. 6.8(e). The damper motion (Fig. 6.8(d))also indicates the motion 
is tending toward a stable mode of response. 

For a rotor speed of 37,000 RPM the response for five cycles is 
given by Fig. 6.9(a) -(f). The rotor still shows a very small amplitude 
in comparison to the journal relative motion which is indeed unstable 
and spirals out in cycles two through five at an average whirl of 0.50. 

The damper motion Is also observed to whirl at large amplitudes as compared 
to the motion at 29,000 RPM. 

If the bush clearance is reduced to 10 mils the stability map indicates 
the system should be unstable. Figs. 6.10(a) -(f) give eight cycles of 
motion and indicate that the journal relative motion has a sustained 
circular whirl of 0.50 while the rotor motion is very small' and the 
damper motion has non-circular motion but similar to the stable case of 
Fig. 6.8. It must be remembered that the stability chart indicates the 
stability at "a point". As soon as the bush motion begins the equilibrium 
position changes and so does the stability of the system for cases near 
the threshold speed. However the journal relative motion in Fig. 6.10(e) 
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FIG. 6.8(a) ABSOLUTE ROTOR MOTION FOR TWO CYCLES (N=29, 000RPM, Ks“500LB/MIL) 

(b) JOURNAL RELATIVE MOTION BELOW STABILITY THRESHOLD 

(c) DAMPER SUPPORT TRANSIENT RESPONSE (CLEARANCE®! 7 , 5MILS ) 
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FIG. 6.9(a) ABSOLUTE ROTOR MOTION FOR TWO CYCLES (N=37, 000RPM, Kq=500,000) 

(b) JOURNAL RELATIVE MOTION ABOVE STABILITY THRESHOLD SHOWING HALF 
FREQUENCY WHIRL 

(c) DAMPER SUPPORT TRANSIENT RESPONSE ( CLEARANCES 7. 5MILS) 






203 


FIG. 6.9(d) ABSOLUTE ROTOR MOTION FOR CYCLES TWO THROUGH FIVE <»"37 t 000RPM, 
Kc:“500LB/MIL) 

(e) ^JOURNAL RELATIVE MOTION SHOWING UNSTABLE HALF- FREQUENCY WHIRL 

(f) DAMPER SUPPORT SUSTAINED WHIRL (CLEARANCE^! 7 . 5MILS ) 
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FIG. 6.10(a) ABSOLUTE ROTOR TRANSIENT FOR FOUR CYCLES (N=29, 000RPM, Ko^OOLB/MIL, 

(b) JOURNAL RELATIVE MOTION IN UNSTABLE REGION OF STABILITY MAP 

(c) DAMPER SUPPORT RESPONSE ( CLEARANCE 3 1 OMILS ) 






2 

c* 

os 

o 

o 

o 

* 

o% 

fM 

« 


o 


3 

S' 


o 

I— t 

w 

X 

a 

o 


« 

k< 

*-4 


B 

Z,'~^ 
M CO 
<! J 

H M 

co £ 

P O 

o w 
zo 

<i 

Sd 

a p 
z o 

Mw 

gw 

O CO 

gp 

O th 
X CO 

w 

g* 

MH 
&* 

K»4W(3k 

HftJP 

im^* 

p sac as 

PPZ w 

Q O Od Ou 

wops: 

STSS 

>? 


OS 

rs 


-3 

8 

0 

1 

t *4 

B 

w 

i 

CS 

H 


0 ) <« 


vC 

o 

M 

w 


205 





does have a sustained whirl even though the limit cycle is very small, 
and verifies the predicted instability. 

When the bush clearance is reduced to 5 mils the same type limit 
cycle motion for the journal is indicated in Fig. 6.11(b) and the bush 
motion of 6.11(c) indicates less stability than for the 10 mil damper. 

To demonstrate how the damper has improved the rotor motion a rigid 
support case for five cycles of motion was computed (see Fig. 6.12). 

The rotor and journal spiral outward in a very unstable whirl orbit after 
the initial drop of the rotor and lift of the journal to achieve the 
proper rotor-journal force balance. 

The flexible rotor (6/c = 1.6) has the same upper threshold as the 
rigid rotor for the larger clearance values. The same speeds of 29,000 
and 37,000 RPM were run for four cycles beginning from the end conditions 
of the rigid rotor case (plus a static bow for the flexible shaft). 

These orbits appear in Fig. 6. 13(d) -(f) where! the journal relative 
motion is observed to go into the half-frequency whirl above the predicted 
threshold speed. The bush motion is considerably larger than in the rigid 
rotor case but the rotor is noted to be seeking another position than the 
one given as a starting condition. Fig. 6.14(a) and 6.14(b) indicate the 
unstable journal motion for rigid supports and the fact that the rotor 
motion has not gone unstable during these five cycles of motion. The 
damper support has improved the response a great deal as compared to the 
whirl orbit obtained on the rigid supports. 

In order that the lower region of instability could be verified 
the 17.5 mil bush damper was tested. The system specifications appear 
in Table 6.1 and the information in Table 6.2 indicates an unstable whirl 
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(b) JOURNAL TRANSIENT SHOWING HIGHLY UNSTABLE SYSTEM 
AFTER INITIAL RESPONSE TO PERTUBATION 
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FIG. 6.13(a) ABSOLUTE ROTOR RESPONSE FOR CYCLES TWO THROUGH SIX (N=29, G00RPM, 
Kq~5LB/MIL) 

(b) 5 JOURNAL RELATIVE MOTION FOR STABLE RESPONSE PORTION OF STABILITY 

MAP c 

(c) DAMPER SUPPORT TRANSIENT (CLEARANCE 8 *! 7 . 5MILS ) 
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FIG. 6.13(d) ABSOLUTE ROTOR RESPONSE FOR CYCLES TWO THROUGH SIX (N=37, 000RPM, 
K S =5LB/MIL) 

(e) JOURNAL RELATIVE MOTION INDICATING AN UNSTABLE HALF- FREQUENCY WHIRL 

(f) DAMPER SUPPORT TRANSIENT RESPONSE (CLEARANCE 0 ! 7 . 5MILS ) 
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FIG. 6.14(a) ABSOLUTE ROTOR RESPONSE (N=29, OOORPM, K S =5LB/MIL, 

RIGID DAMPER SUPPORT) 

(b) JOURNAL TRANSIENT SHOWING UNSTABLE HALF- FREQUENCY WHIRL 





STABILITY ANALYSIS OF THE FLOATING BUSH DAMPER 

ROTOR SPECIFICATIONS 

M2 - ** 0.0 0 LB. CS * 0.00 LB-SEC/IN 3 * 0.0 LB/IN 

<S V 5000* 0 LB/IN. IC "* Q.OO' LB -SEC/ IN 

ORDER = 12 
OCR * 12 


JOURNAL BEARING SPECIFICATIONS 

'ij i 5. 0 Of LBV RJ S i. Q I N . MUJ =" 1.0 QBE- 05 REYNS" 

CJ *_ .0050 IN. L J ■». 1.0 IN. 


FLOATING BUSH SPECIFICATIONS. ' 

MB s 2.50 LB. RD = 1.20 IN. MUS * l.OOOE-05 REYNS 

CB = i3175 "IN. LB ii “l.Oa IN. 

PRELOAD-X s 0.00 
PRELOAD-Y * 0.00" 

K1XX = 5772.0 LB/IN 

"K1YY = 5772.0 LB/IN ' ' ““ 


STEADY-STATE POSITION 


ROTOR' 
X2S = 

.743* MILS 

Y2S * 

-12.5096 HlLS 

DLTA = 

1.6000 

JOURNAL 

XRJ = 

.7436 NILS 

YRJ = 

- . Pi 5F“RTU5 

BUSH 

■XB.S ‘ * 

D • 00 GO ruLS 

YBS = 

-4.7644 MILS 


.0175 


3.500 


A 7 74 



WNTBI — KCUNT 
0 0 


TABLE 6.1 FLOATING BUSH DAMPER ROTOR SYSTEM SPECIFICATIONS 
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of 0.35 when the running speed Is 4775 RPM. This information was obtained 
from the stability program SDSTB. The transient program response plots 
are presented in Fig. 6.15(a) - (!) where cycles 0-2 were started from 
approximate steady-state conditions and show no large transients. Figs. 
6.15(d), (e), (f) for cycles 2-6 indicate a large whirl developing at 
the rotor and support with the journal loosing some of its load capacity 
(Fig. 6.l5(e))due to the whirl of the bush as indicated in Fig. 6.15(f). 
Cycles 6 through 10 indicate a sustained whirl of the rotor and bush at 
slightly over one-third as shown by the timing marks (this was as pre- 
dicted by SDSTB - Table 6.2). Note that the journal relative motion is 
actually very "stable" In this mode of operation (Fig. 6. 15(h)) and is 
thus tracking out the same motion as the bush damper ie., ~ 1/3 whirl 
rate ) . 

For a speed ratio of 4.52 ( 1 2,000 RPM) the system is stable as 
indicated by Fig. 6.7. The transient response for both a large initial 
displacement of the journal and a very small pertubation were calculated 
and are given in Figs. 6.16 and 6.17. For a large pertubation the damper 
response is not readily reducing. The journal decreases initially 
(Fig. 6. 16(b)) but has a sustained limit cycle due to the bush transient 
motion. Fig. 6.16 (k) indicates that after 16 cycles the journal relative 
motion is very slowly reducing in amplitude whereas the rotor and bush 
(support) motion has a relatively large motion. However these are not 
in phase as they were for the case of 4775 RPM (unstable condition). 

Figs. 6.17(a)- (f) indicate eight cycles of motion at 12,000 RPM for 
a very small pertubation. The motion is considerably more stable than 
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FIG. 6.15(a) ABSOLUTE ROTOR MOTION FOR TWO CYCLES (N=4,775RPM, Kg^LB/MIL) 

(b) JOURNAL TRANSIENT 

(c) DAMPER SUPPORT TRANSIENT (CLEARANCE - *! 7 . 5M1LS ) 
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FIG. 6.15(d) ABSOLUTE ROTOR TRANSIENT WHIRL FOR CYCLES TWO THROUGH SIX (N- 
4,775RPM, Kq=5LB/MIL) 

(e) JOURNAL RELATIVE TRANSIENT SHOWING NEW RELATIVE EQUILIBRUM POSITION 

(f) DAMPER SUPPORT TRANSIENT MOTION (CLEARANCES 7 . 5MILS ) 
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FIG. 6.16(a) ABSOLUTE ROTOR MOTION FOR FOUR CYCLES (N=l 2 , 000 RPM, Ko«5LB/MIL) 

(b) JOURNAL RELATIVE MOTION FOR LARGE INITIAL DISPLACEMENT PERTUBATION 

(c) DAMPER TRANSIENT RESPONSE ( CLEARANCE® 1 7 . 5MILS ) 
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12,000, K S =5LB/MIL) 

JOURNAL SUSTAINED LIMIT CYCLE 
DAMPER TRANSIENT WHIRL 
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FIG, 6,16(j) ABSOLUTE ROTOR MOTION FOR CYCLES TWELVE THROUGH SIXTEEN (N= 
1 2 , OOORPM, Ko=5LB /MIL ) 

(k) JOURNAL ORBIT 

(l) DAMPER SUSTAINED WHIRL TRANSIENT MOTION 
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FIG. 6.17(a) ABSOLUTE ROTOR MOTION FOR FOUR CYCLES (N=12,000, K s »5LB/MIL) 

(b) JOURNAL TRANSIENT FOR SMALL INITIAL PERTUBATION 

(c) DAMPER SUPPORT TRANSIENT ( CLEARANCES 7 , 5MILS ) 
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for the large per+uba+ion but the bush damper and rotor are giving some 
Initial motion. The motion would eventually reduce to the steady-state 
equilibrium position as predicted by the stability analysis. 

The sustained transient due to the large pertubation in what should 
be a stable system is the greatest drawback to squeeze damper design. 

This type of performance could be catastrophic in an actual rotor system. 
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CHAPTER V I I 


CONCLUSIONS 

7.1 The Importance of Analytic Simulation 

Due to the ease of manufacture and Installation of plain journal 
bearings in pumps, transmissions, and small motors, there is need for 
the dispersion of information on the stability characteristics and 
unbalance response of plain journal bearings. Plain journal bearings 
operating under reduced-load conditions can exhibit large whirling and 
eventual failure of the bearing surface caused by the cyclic loading 
from the whirling journal. Isolated journal bearing performance studies 
can produce this much needed stability data and unbalance response in- 
format i on . 

Manufacturers of turbine and compressors conduct test runs on actual 
equipment to verify the design and performance of their products. Each 
test unit may represent from $200,000 to In excess of a half-million 
dollars in hardware and instrumentation. A single bearing failure or 
rotor-shaft instability could destroy the entire test facility in a matter 
of seconds if the unit malfunction were not detected or the worst case 
would be if the test group could not recognize the presence of the in- 
stability and its eventual growth leading to a system failure. 

The instrumentation necessary to observe the rotor shaft and housing 
motion is readily available and easily installed in existing test rigs. 

The complete shaft orbit should be observed to prevent the possibility 
of choosing the incorrect plane to observe the greatest motion. Systems 
having asymmetric bearing characteristics may exhibit instability in one 
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plane and stability in a plane at right angles to first. The incorporation 
of timing marks on the rotor orbits by the use of a reference mark on the 
shaft and a keyphasor probe is also of importance in the determination of 
the rate of whirl. Subharmonic and superharmonic whirl is easily dis- 
tinguished if timing marks are placed on the whirl orbit trace. 

The analytic simulation is important due to the fact that controlled 
excitation can be introduced in the model and the reaction observed from 
that isolated excitation. Combinations of external forces, unbalance, and 
internal damping can also be studied under controlled conditions and 
their response characteristics observed. In an actual system the exact 
form of excitation might be very hard to determine and the only informa- 
tion available is +he response of the unit. The machine analyst must be 
able to infer from the performance of the machine the likely forms of ex- 
citation and the approximate magnitude. The easiest way to obtain this 
recognition ability is by the study of the controlled digital simulation 
of accurately simplified rotor-bearing models. 

7.2 Summary of Major Results and Conclusions 

This analysis has developed and presented results of steady-state 
response, transient simulation, and stability of rigid and flexible rotor 
models. The transient models developed included both linear and nonlinear 
bearing and support characteristics. The results and conclusions may 
be summarized as follows: 

(I) Nonlinear transient analysis of complex rotor-bearing systems 
can be readily produced by digital computation by the approach 
presented in Chapter II. 
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(2) Simplified steady-state models of rotor systems have been 
presented and may be used in the design and analysis of rotating 
equipment in addition to providing initial conditions for the 
transient computer codes. 

(3) Simplified rotor model transient computer codes have been de- 
veloped to check the major transient program of this analysis 
(MODELJ) and are themselves a valuable tool to the machine 
analyst due to the simplicity of the input specifications. 

(4) Transient simulation Is necessary to verify the rate of growth 
of instabilities in nonlinear systems as predicted by linearized 
stabi I ity studies. 

(5) The stability boundaries as given by Choudhury and Gunter (69) 
for the elastic mounted plain journal bearing have been verified 
by transient simulation. 

(6) Unbalance was shown to be advantageous in vertical plain journal 
bearings having elastic support structures. This was previously 
reported for the rigid support vertical journal bearing. 

(7) Phase angle measurements may be used in addition to amplitude 
response to analyze machine performance. The timing mark on 
elliptical orbits can be related to circular synchronous motion 
by graphical construction. 

(8) It has been demonstrated that it is feasible to use the direct 
expansion approach on a 6 x 6 quadratic determinant to deter- 
mine the stability of several important rotor models which had 
not been reported in the literature by this approach to date. 
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The key to time saving in the expansion is a zero check to 
avoid many needless operations. 

(9) The phenomenon of passing through regions of instability has 
been demonstrated by the stability analysis of damper bearings 
and verified by transient simulation. 

(10) It was shown that floating bush dampers indicated as stable in 
linearized stability criterions can have large sustained motion 
when large perturbations are imposed on the system. 

(11) Numerous example rotors have been investigated and illustrate 
the versatility of the simulation approach to design and 
diagnosis of machine failures. 

(12) Damped flexible supports can be used to increase the stable 
operating speed range of high-speed rotors. A properly de- 
signed support reduces the System transmissibi I ity and maintains 
low level response for unbalance excitation in the stable 
operating speed range. 

(13) Suppression of complex aerodynamic excitation can be studied 
for elastic mounted rotor-bearing systems by the stability 
program developed for the floating bush damper bearing. 

(14) Timing marks on simulation orbits correspond to keyphasor 
timing marks placed on actual machines and are excellent 
means of helping to interpret the resulting orbits. 

(15) The more complex transient programs are more useful for 
analysis of existing machine designs and evaluation of proposed 
system modifications than for producing extensive general 
design criteria. 
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(16) Steady-state response i n format! on will be essential for large 
complex rotor simulations to avoid costly undesired transient 


behavior. Zero initial conditions may be acceptable for impact 
unbalance studies but instabilities due to aerodynamic excita- 
tion, internal damping, and fluid-film bearings are best studied 
from steady-state initial conditions with small pertubations. 

7.3 Suggestions for Future Research 

This analysis has developed several computer codes for the simulation 
of rotor-bearing systems. Additional design criterion can be produced 
from the rotor models presented and will be of extreme interest to turbine 
and compressor manufacturers. Some of the areas that need extension are 
as fol lows: 

(1) Extensive studies of rotor acceleration rates on forces trans- 
mitted and transient response. 

(2) Inclusion of skewed disk excitation on the flexible rotor 
simulation for studies of transient behavior. 

(3) Analysis of coupled lateral -torsional modes of vibration for 
the flexible rotor, 

(4) Extensive stability maps produced showing the effect of damper 
supports on aerodynamic and internal friction instability. 

(5) Stability and transient response of the flexible rotor-bearing 
system on closed-end pressurized damper supports. 

(6) Verification of the response characteristics of the damped 
support by controlled experimental models of flexible rotor- 
bearing systems. 
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( 7 ) 


Verification of the reduction of the transmissi bi I ity factors 
by experimental techniques. 

(8) Extensive transient orbit studies of the response of the non- 
linear fluid-film dampers to large perturbations. 
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